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I.  INTRODUCTION. 

This  report  is  a detailed  analysis  of  the  mechanical  interaction  between  a 
clip  strain  gage  and  a solid  propellant  uniaxial  tensile  sample  and  the 
effects  of  this  interaction  on  the  strain  measurement  provided  by  the  gage. 

The  analysis  was  first  undertaken  in  the  summer  of  1976  to  assess  an  apparent 
potential  technical  problem  which  was  then  developing  on  AFRPL  Contract 
FO^ 1 1 -76-C-OOAl , "Improved  Solid  Propellant  Mechanical  Properties  Measurement 
for  Structural  Analysis  Input"  (hereinafter  referred  to  as  the  Improved 
Properties  program),  with  Thiokol  Corporation/Huntsville  Division.  The 
planned  testing  effort  under  that  program  included  high  rate  tensile  tests 
with  a closed,  liquid-filled  pressurization  chamber  enclosing  the  sample. 

A key  mechanical  response  parameter  to  be  determined  from  these  tests  was 
Poisson's  ratio;  the  contractor  intended  to  determine  this  property,  under 
dynamic  pressurization  and  straining  of  the  sample,  from  separate  measurements 
of  axial  and  transverse  strain  in  the  test  specimens.  These  strains  were 
to  be  measured  using  clip  strain  gages. 

Believing  the  fabrication  of  highly  compliant  gages  to  be  a relatively 
straightforward  problem,  the  contractor  had  planned  no  detailed  analysis  of 
the  measurement  problem  to  define  acceptable  values  for  the  gage  stiffness. 

As  the  Air  Force  Project  Manager,  the  author  decided  to  independently  analyze 
the  problem  so  that  criteria  would  be  available  if  needed.  The  results  of  this 
rudimentary  analysis,  together  with  laboratory  tests  of  gaged  propellant  samples 
conducted  by  the  contractor,  led  to  the  mutual  conclusion  that  the  high  pre- 
cision required  for  Poisson's  ratio  determination  could  not  be  achieved 
(within  the  resources  available)  with  clip  strain  cages.  The  approach  was 
abandoned  at  that  point.  (The  contractor  eventually  devised  an  acceptable 
alternate  technique  for  Poisson's  ratio  determination,  which  is  discussed  in 
Reference  1.) 
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The  early  analysis  results  agreed  qualitatively  with  the  gaged  sample  test 
results,  but  substantially  over-predicted  the  measurement  bias.  Some  of  the 
assumptions  used  in  the  analysis  were  re-examined  and  found  to  be  in  error. 
The  analysis  was  corrected  and  close  agreement  was  obtained  with  the  test 
results.  A comparison  with  check  problems  using  the  TEXGAP-3D  finite  element 
code  also  confirmed  the  revised  assumptions. 


Since  the  results  of  this  analysis  were  judged  to  be  of  value  to  the  solid 
rocket  propulsion  community  at  large,  the  decision  was  made  to  publish  a report. 
In  the  process  of  writing  the  report,  many  areas  of  further  analysis  were 
pursued,  both  to  strengthen  the  basic  structure  of  the  analysis  and  to  explore 
the  overall  propellant  strain  measurement  problem  in  greater  detail. 


2. 


SUMMARY. 


The  analyses  described  in  this  report  provide  an  overall  model  for  the  mechani- 
cal interaction  between  a clip  strain  gage  and  a uniaxial  propellant  sample 
of  a particular  configuration.  The  interaction  analysis  separately  considers 
the  effects  due  to  "beam  deformat  ion"  of  the  propellant  sample  under  the 
clip  gage  resisting  force  and  the  local  deformation  of  the  gage  mounting  pins 
under  the  same  force.  These  two  components  of  the  interaction  are  then 
combined  to  determine  the  net  interaction  effects  as  a function  of  gage 
st i f fness. 

Typical  clip  gage  conf igurations  are  analyzed  in  a parametric  fashion,  and 
the  effects  of  gage  design  factors  (geometry  and  materials)  on  gage  output 
and  (through  the  gage  stiffness)  on  relative  strain  measurement  bias.  The 
analysis  shows  that  the  original  Thiokol  prototype  gage  produces  a very 
large  strain  bias  (roughly  16  percent  of  the  applied  strain  at  a propellant 
modulus  of  500  psi)  and  that  the  final  Thiokol  prototype  gage  shows  a sub- 
stantial improvement  (to  about  6 percent  of  the  applied  strain  at  500  psi 
propellant  modulus).  Since  this  level  of  measurement  bias  is  not  desirable, 
approaches  to  reducing  the  measurement  bias  are  discussed  and  analyzed.  It 
is  shown  that  major  reductions  in  bias  (without  sacrificing  gage  output)  are 
potentially  feasible. 

The  elastic  analyses  used  to  accomplish  the  above  results  is  then  extended  to 
include  viscoelastic  propellant  behavior,  obtaining  the  time  dependent  clip  gage 
response  to  a series  of  step  strain  increments  applied  to  the  propellant  sample. 
These  results  are  shown  to  compare  extremely  well  with  actual  measurements  made 
by  Thiokol  on  a gaged  propellant  sample.  A limited  study  of  gage  response 
under  strain  cycling  and  short  time  effects  of  clip  gage  prestrain  is  also 
presented. 

The  final  conclusions  are  that  clip  gage  improvements  are  needed  to  reduce 
measurement  bias  due  to  gage-propellant  interaction,  that  such  improvements 
are  possible,  and  that  the  interaction  problem  can  be  analyzed  with  acceptable 
accuracy . 
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3.  ANALYSIS  OF  MECHANICAL  CONTRIBUTIONS  TO  MEASUREMENT  ERROR. 

This  section  comprises  the  main  body  of  the  report.  The  major  contributions 
to  measurement  error  produced  by  the  mechanical  behavior  of  a system  con- 
sisting of  a uniaxial  propellant  sample  , a cl  ip  strain  gage,  and  the  mounting 
pins  which  connect  the  clip  gage  to  the  sample  are  analyzed  in  detail.  The 
analyses  reported  here  developed  in  a rather  unplanned  fashion.  It  was 
evident  that  the  first  simple  analyses  (section  3.1,  which  treats  the  pro- 
pellant sample  as  a beam  with  a superimposed  axial  load,  and  3.**,  which 
consists  of  analyses  of  the  clip  strain  gage  structural  behavior  and  its 
interaction  with  the  sample  as  characterized  in  3.1)  ignored  a very  important 
souce  of  error.  This  source  of  error,  local  deflection  of  the  mounting  pins, 
was  therefore  analyzed  separately  in  3-2  and  the  combined  effects  (beam 
deformation  plus  local  deflection)  were  obtained  as  described  in  3*3.  Finally, 
the  elastic  analyses  of  3-1,  3-2,  and  3-3  were  modified  to  incorporate 
time-dependent  propellant  modulus  to  produce  the  results  discussed  in  3*5. 

3. 1 Gage-Sample  Interaction  Due  to  Beam  Deformation  of  Sample. 

The  analyses  in  this  subsection  were  conducted  to  estimate  the  direct  effects 
of  gage  stiffness  on  the  bias  of  strain  and  force  measurements  on  a gaged 
uniaxial  tensile  propellant  sample.  The  sample  configuration  used  in  the 
analysis  is  an  end-bonded  sample,  4 inches  long,  with  a I/2-inch  by  1/2- 
inch  square  cross  section.  This  sample  configuration  is  the  one  used  for 
high  rate  uniaxial  tests  in  the  Improved  Properties  program. 

The  analyses  which  follow  are  based  on  simple  elastic  beam  theory  and 
include  all  the  assumptions  thereof  (plane  sections  remain  plane,  etc.). 

By  far  the  most  important  assumption,  of  course,  is  linear  elasticity. 

While  propellant  is  not  a linear  elastic  material,  it  is  often  successfully 


1.  In  this  context,  zero  bias  in  attained  when  the  stress  and  strain  dis- 
tributions in  the  gaged  sample  are  identical  to  those  which  would  occur  if 
the  gage  were  absent. 
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The  structural  model  analyzed  in  this  subsection  is  shown  in  Figure  1. 

While  the  ends  of  the  propellant  specimen  are  not  truly  fixed,  the  fixed 
ends  of  the  model  more  nearly  simulate  the  true  end  constraints  than  do 
pinned  ends.  The  effect  of  the  loads  P is  statically  equivalent  to  that 
of  two  loads  of  the  same  magnitude  acting  at  the  sample  centerline  (axial 
load)  plus  a pair  of  moments  of  opposing  directions  of  magnitude  = Ph/2. 
(The  effect  of  load  application  at  some  height  e above  the  surface,  as  shown 
in  Figure  2,  is  to  increase  the  moments  by  an  amount  AMq  = Pe.)  These 
forces  and  moments  are  assumed  to  be  uniformly  distributed  over  transverse 
sections  of  the  beam  which  include  the  points  of  application  of  the  loads  P, 
so  that  simple  beam  theory  can  be  used.  This  assumption  introduces  an 
unknown  (but  hopefully  small)  error,  since  the  gage/sample  interactions 
involve  local  deformations  near  the  mounting  surface  as  well  as  the  aver- 
age deformations  of  the  beam  cross  sections.  The  local  deformations  be- 
come less  important  as  the  area  of  the  clip  gage  mounting  pads  (Figure  2) 
is  increased.  (The  local  deflections  for  a different  mounting  system  using 
straight  pins  piercing  the  sample  are  analyzed  in  3-2). 
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treated  as  such  by  using  “quasi-viscoelastic"  analysis  procedures.  For  the 
purposes  of  the  present  investigation,  the  use  of  elastic  theory  appears  to 
be  justified.  The  elastic  modulus  values  used  in  the  analyses  should  be 
considered  as  effective  instantaneous  values  of  either  the  relaxation  modulus 
or  the  constant  strain  rate  modulus,  depending  on  the  test  being  considered; 
both  types  of  loading  are  planned  for  these  samples.  It  Is  obvious  that  a 
linear  elastic  analysis  cannot  account  for  any  effects  of  modulus  dependence 
on  strain;  therefore,  since  most  propellants  show  strain  dependence,  the  use 
of  the  analysis  results  to  correct  for  measurement  bias  would  be  somewhat 
questionable.  However,  the  analysis  should  become  more  accurate  when  the 
perturbations  to  the  "ungaged"  sample  behavior  become  small;  therefore  the 
analysis  results  can  be  reasonably  applied  to  specify  conditions  for  con- 
trolling bias  to  low  levels. 
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Certain  known  mechanical  sources  of  measurement  bias  were  not  analyzed. 
These  are  listed  here  to  prevent  any  appearance  that  all  sources  of  direct 
mechanical  error  are  covered  by  the  analysis,  even  within  the  context  of 
linear  elastic  analysis: 


Nonuniform  axial  stress/stra in  distribution  in  the  sample 
under  uniaxial  load  (in  the  absence  of  clip  strain  gages) 
produced  by  end  constraints.  (See  discussion  in  3 - 5 • ) 


b.  Additional  perturbations  caused  by  other  clip  strain  gages 
(e.g.,  transverse-strain  measuring  gage)  which  may  be  present 
in  addition  to  the  gage  included  in  the  model. 

c.  Deviations  from  the  assumed  distribution  of  the  clip  gage 
restraining  force  across  the  sample  surface  (except  for  the 
straight  pin  mounting  system  analyzed  in  3-2). 

d.  The  additional  restraint  on  the  propellant  surface  deforma- 
tion produced  by  the  mounting  pads  (except  for  the  straight 
pin  mounting  system  analyzed  in  3-2). 

All  of  these  factors  produce  measurement  bias  which  adds  to  the  bias  con- 
sidered in  this  analysis.  Every  effort  should  be  made  to  minimize  their 
effects  in  designing  a clip  strain  gage  and  its  mounting  system. 


The  analyses  on  the  following  pages  explore  the  effects  of  clip  strain 
gage/sample  interaction  due  tc  beam  deformation  on  the  bias  in  strain 
measurements  using  the  gage  itself,  in  strain  measurements  based  on 
optical  measurement  of  the  motion  of  reference  points  on  the  sample,  and 
in  stresses  determined  from  direct  force  measurements. 


1 


3*  !• I Bias  in  Strain  Gage  Measurement. 


The  effects  of  the  axial  load  and  the  moment  will  be  considered  separately 
and  then  combined.  For  the  axial  load,  the  stress2  between  z - a and  z - b 
in  Figure  2 is  (Reference  I,  p.  82): 


(Ao  ) 

z axial 


and  the  strain  Is 


(Aez)  axial  m I (Aoz  + vAox  + vAoJ 


z x y axial 


(Acz)  axial 


since  Ao  » Ao  » 0. 

^ y 


For  the  applied  moments^,  the  reaction  force  at  z « 0 is  (Reference  2,  p.  113): 


R,  * ~ p Hat  - a2)  - (be  - b2)] 


R,  » -p1  [(ae  - a2)  - (be  - b2)]. 


2.  It  is  understood  that  the  sample  is  under  some  yet-to-be-determined 
stress  and  strain  as  a result  of  the  applied  deformation  at  the  sample  ends 
and  that  these  stress  and  strain  values  are  additive  components  produced  by 


P and  M . 
o 


3.  Note  that  beam-column  Interactions  are  ignored  in  the  ensuing  dlsccusion. 
The  impact  of  this  omission  has  not  been  examined. 
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However,  substitution  of  b - i-a  produces  the  result 


R,  *0 


(3) 


Also,  the  reaction  moment  at  z - 0 is  (Reference  2,  p.  113): 

Ml  “ ~T  t(Ha-3a2-t2)  - (l*tb-3b2-t2)  ] 

j r 

or,  substituting  Mq  * Ph/2  and  b “ t-a, 


M.  - ^ Uta-S.2]. 
1 2i 


(4) 


The  bending  moment  between  z * a and  z * b is  (Reference  2,  p.  113): 


M = -M.  + R.z  + M 
i o 


(5) 


and  the  resulting  flexural  strain  between  z « a and  z * b is 


(Ac  ) . » ■=■  (A o + vAa  + vAa  ) . 

' z bending  E z x y bending 


so 


Uez) 


bending 


bending 

E 


M(h/2) 

El 


(6) 


The  actual  quantity  of  interest  is  the  average  strain  over  the  strain  gage 
measurement  base  (gage  length),  i.e.  between  z ■ a and  z * b.  However, 
since  Rj  is  zero  in  Equation  (5),  the  strain  is  constant  in  this  region  and 
the  average  is  equal  to  this  constant  value.  A similar  observation  applies 
to  the  axial  strain  component  given  by  Equation  (2). 


For  the  sample  dimensions  given  in  Figure  1,  the  cross  section  moment  of 
inertia  is 


I 


1 

12(2)4 


1 

192 


\ 


and,  with  h « 1/2,  Equation  (6)  becomes 


“'/bending  ‘ 48  H/E' 


(6a) 


For  a clip  strain  gage  length  of  Lq  “ O.^fOO  inches  and  the  sample  dimensions 
of  Figure  1,  Equation  (2)  yields 


(Aez)axial 


3.60 


E 


and  Equation  (*»)  yields 


Mj  - +0.0250  P. 

Therefore,  from  Equation  (5), 

M - -0.0250  P + 0 + - 0.2250  P 

since  h - 0.5. 

Substitution  of  this  value  of  M into  Equation  (6a)  results  in 
(Aez) bending  = A8  (°’225°)  p/E  = ,0*80  p/E. 

Final ly, 

Ac  - (Ae  ) . . + (Ae  ),  , . 

z z'axial  z bending 

= 3.60  P/E  + 10.80  P/E  - 1 40  P/E  (7) 

or  (Ae  )/e  - (14. 40  P )/(e  E)  • strain  attenuation  /g\ 

x o o o gage  strain  ' ' 

where  e is  a specific  value  of  gage  strain  (AL  /L  ) 
o o o 


and  Pq  is  the  reaction  force  at  each  gage  foot  corresponding  to  the  strain 
cq.  Since  the  local  sample  strain  without  the  gage  would  be  cq  + Aez>  the 
relative  measurement  bias  for  strain  is 


-Ae 


-Ae 


'Re  sample  strain  e + Ae 

o z 


or,  from  Equation  (8), 


■14.40  P 


-14.40  P /E  ...... 

o o 

eRe  " e + 14.40  P 7e  “ El  + 14.40  P 
o o o o 


(9a) 


For  e ■ 0.040  in/0.400  in  » 0.1, 
o 

„ . -14.40  P 

gage  strain  - sample  strain  = o 

eRe  sample  strain  0.1  E + 14.40  P 


(9b) 


Calculations  were  carried  out  for  different  combinations  of  P and  E for 

o 

Lq  = 0.400.  The  results  are  shown  in  Figure  3. 


The  above  calculations  were  carried  out  for  a gage  length  of  0.400  inches. 

To  examine  the  effects  of  changes  in  gage  length  on  the  strain  measurement 
bias,  Equation  (9)  can  be  solved  for  P , yielding 

P ~eRe  (Eeo)  (10) 

o 14.40  (I  + eD  ) 

Re 


4.  Note  that  if  the  gage  preload  (force  at  zero  sample  strain)  is  zero, 
PQ/eo  for  any  gage  strain  eQ  is  a constant  for  the  gage.  Zero  preload  was 

assumed  for  the  analyses  in  this  subsection  and  in  3-3.  If  a preload 
exists,  P,-/e0  becomes  a function  of  applied  strain.  This  added  complication 

is  addressed  in  3-5. 
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Ejl.:i(S.tjt'iijrn  Gagti)  itrtalh  Measurement  for 


Gaged  Tens i To  Spice i man; 


!*i  Loriytfi 


and  K can  be  substituted  for  14.40,  where  K is  dependent  on  lQ  and  M in 

accordance  with  Equations  (2)  and  (6).  K can  then  be  re-evaluated  in 

terms  of  different  values  of  L and  different  values  of  loading  eccentricity 

o 

(i.e.,  e i<  0 in  Figure  2).  These  calculations  were  performed  for  PQ  at 
-1  percent  bias  and  the  results  are  shown  in  Figure  4 (effect  of  gage  length 
variation)  and  Figure  5 (effect  of  eccentricity  variation).  In  drawing 
conclusions  from  Figure  4,  the  effect  of  changes  in  gage  length  on  the 
actual  value  of  Pq  for  the  clip  strain  gage  must  be  considered. 


3.1.2  Bias  in  Optical  Strain  Measurements. 

An  alternate  means  of  measuring  the  sample  strain  under  appropriate 
laboratory  conditions  is  to  visually  observe  (e.g.  with  a cathetometer) 
the  displacement  of  a pair  of  reference  marks  on  the  propellant  surface. 
The  result  is  a measurement  of  the  average  strain  over  the  visual  gage 
length  L 1 (see  Figure  1). 

The  axiai  components  of  stress  and  strain  are  given,  respectively,  by 
Equations  (i)  and  (2)  for  the  portion  of  the  sample  between  z * a and 
z - b (Figure  1).  For  values  of  z less  than  z = a,  it  can  be  shown  that 


^°z^axial , 0 <_  z <_  a A ( *,)' 
^Aez^axial , 0 <_  z <_  a AE  ( £/” 


The  average  axial  strain  over  the  gage  length  L 1 is  then 


(Aez)  axial 


f.  * - fr  fh  ft)  - fw  (4) 

° J /*- L '\  ° J /l-L  *\  J 

i-^j  L M * 
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RROOUCE 


GAGE 


Fiigurti 


'RdvjiJ.i  red 


<iii  Force 


I'd  i n Gni:e) 


itQ.Pfodiice:';a 


■Maavyrod;  St 


&j»r,i$d’  Tcfti 


rmen 


T 


(Ae  ) . . 

z axial 


2P  | 

L ' AE  I 1 


Ji-L. 


i-i 


2(M 


(13) 


It  can  be  demonstrated  that  Equation  (13)  reduces  to  Equation  (12)  for 


L . 
o 


The  bending  moment  for  values  of  z between  z ■ a and  z = b is  given  by 
Equation  (5).  For  values  of  z less  than  z = a,  the  bending  moment  is 
(Reference  2,  p.  113): 


M-  -Mj  +R)x=-M),  0 £ z < a. 

The  average  flexural  strain  is  therefore 

1/2 


(Ae  J 


z bending  L 


e/2 

Mdz 


m & 


2c  I /*Ph  .2v 

' Lj'EI  -2  > d2  * 


.e/2 

Ph 


2e 


2 (2  ea- C2 ) + p- 


dz 


C Ph 
E I L 1 „2 

o e 


e (e-L  -) 

+ eaL  - e a 

o 


T 


(14) 


For  l 


Lo,  Equation  1 A reduces  to  Equation  (A)  if  c ■ h/2. 


Substitution  of  dimensions  from  Figure  (I),  with  L ' = 3 inches,  into 
Equations  (13)  end  (I  A)  produces  the  results: 


(15) 

lJrz>l»«^ln9  ■ r • 

(16) 

The  most  likely  values  for  C are  +h/2,  0,  and  -h/2  (i.e.,  the  top  surface, 
neutral  surface,  and  bottom  surface  of  the  specimen  as  shown,  respectively, 
in  Figure  I).  Using  these  values  for  C,  the  following  results  are  obtained 
for  the  average  total  surface  strain  over  the  3" inch  gage  length,  where 

Ae  “ I 


) . i + (Ac  ),  : 

axial  z bending 

Top  Surface 

Ac 

z 

*=  1.2667 

P 

E 

(17) 

Neutral  Surface 

Ae 

z 

= 0.8267 

P 

E 

(18) 

Bottom  Surface 

Ae 

z 

= 0.4267 

P 

E * 

(19) 

Now  consider  the  effect  of  a load  Pq  which  corresponds  to  a clip  strain  gage 

strain  e . For  the  specific  (strain  gage)  gage  length  and  sample  dimensions 
o 

used  in  deriving  Equations  (17)  through  (19),  the  strain  attenuation  was 
previously  shown  to  be  given  by  Equation  (7): 


Ae  = I 4. 40  p — 
z E 


so  that  the  sample  strain  without  the  gage  present  would  be 

P 

e + Ae  = e + 14.40  . 

o z o E 

The  strain  measurement  bias  over  the  gage  length  L ' = 3 inches  is  therefore 


e'  » relative  bias 
Re 


-Ae 


Lo'  = 3 incheS  e + 14.40  ^ 


or,  expressing  Equations  (17) » ( 1 8) , and  (19)  in  the  form  Aez  = Kc  and 


substituting  into  the  above  equation, 
K P 


Re 


c E + 14.40  ^ 
o E 


_ gage  length  strain  - sample  strain  ,^q) 

P sample  strain 
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Comparison  of  Equation  (20)  with  Equation  (9)  shows  that 


!Re  TOO 


where 


Kc/14.40  = 0.0880  (Top  Surface) 


= 0.0574  (Neutral  Surface) 

= 0.0296  (Bottom  Surface) 

Thus  the  bias  introduced  in  the  strain  as  measured  by  direct  visual  observa- 
tion of  the  gage  length  L 1 (because  of  the  presence  of  the  clip  strain  cage) 
will  vary  from  2.96  percent  to  8.80  percent  of  the  bias  introduced  into  the 
clip  gage  measurements  themselves  (gage  length  = 0.400  inches).  The  specific 
proportion  depends  on  the  location  chosen  for  the  optical  gage  length  marks 
(top,  bottom,  or  neutral  surface).  It  appears  that  the  best  location  for 
the  marks  would  be  the  bottom  surface  since  the  least  bias  would  result  from 
that  choice. 

3.1.3  Bias  in  Stress  Measurement. 


The  increment  in  measured  stress  caused  by  the  clip  gage  resisting  force  P is 
given  by  Equation  (11): 


v axial , 0 <_  z <_  a A \S,  / 


or,  for  Lq  - 0.400,  i =*  4,  and  A = (1/2)  = 1/4: 


or  (Ac  ) . , n 1 ” 0.400  P 

z axial,  eo  * 0. I o 


Now  from  Equation  (9), 


(sample  strain)  , = (c  ) . 

' r ungaged  z ungaged 


IWp 

o 

oTTTTTOO" 


and,  since  o ■ Ee  , 
z z 


7)  ungaged  = 0.1 


0.1  E 

i 4 . 4o  p 

o 

0.1  E + 14.40  P 


0.1  E (0.1  E + 14.40  P ) 
o_ 

0.1  E 


(Ao  ) , ^ - 0.1  E + 14.40  P^ 

z ungaged  £ = C.l  o 

and,  using  Equations  (22)  and  (23): 


°z  axial , e =0.1 
. . . . ' o 

e„  = relative  stress  bias  = -r—r i 

Ra  rV  ungaged  eo  = 0.1 


so  the  relative  stress  bias  is 

0.400  P 
o 

CRo  “ 0.1  E + 1 4. 40  P 


The  bias  produced  by  different  combinations  cf  PQ  and  E for  Lq  = 0.400  is 
shown  in  Figure  6.  Note  that  this  bias  in  stress  (relative  to  stress  in  an 
ungaged  sample)  is  much  smaller  than  the  bias  in  strain  measurement  using  the. 
clip  gage  (Figure  3)*  The  bias  in  values  of  modulus  determined  from  stress 
and  strain  values  is  therefore  dominated  by  (and  approximately  equal  in 
magnitude  to)  the  strain  measurement  bias.  The  modulus  based  on  uncorrected 
stress  and  strain  measurements  would  be  higher  than  the  true  effective  pro- 
pellant modulus;  for  instance,  if  E * 500  psi  and  Pq  = C.l  lb  (at  eQ  = 0.100), 
the  following  results  are  obtained,  using  Equation  (7)  and  Equation  (24): 
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Measured  Strain  * e - 0.1000  In/In 

P 

Ungaged  Sample  Strain  - eQ  + 14.4 0 -2.-  0.1029  in/ln 

t'n gaged  Sample  Force  ■>  Ee  - (500)  (0. I029)«*5I . 45  lb 

Measured  Force  - 51.45  (1  + eD  ) ■ 51.49  lb 

RO 

Calculated  Modulus  - - 514.9  psl. 

The  bias  In  the  measured  modulus  value  Is  (514. 9-500) /500  - 0.298  or  +2.98 
percent;  inspection  of  Figure  ? shows  a strain  measurement  bias  of  -2.8 
percent  for  the  specified  values  of  E and  Pq,  which  would  produce  a +2.9 
percent  error  in  the  calculated  value  of  E. 
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1 
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3. 2 Gage-Sample  Interaction  Due  to  Local  Pin  Deflection. 

The  analyses  in  this  subsection  were  undertaken  to  estimate  the  added  bias 
caused  by  local  deflections  associated  with  an  alternate  method  of  mounting 
the  clip  strain  gage  to  the  propellant  sample.  In  this  alternate  epproach, 
the  mounting  pad/mounting  pin  system  of  Figure  2 is  replaced  by  straight 
steel  pins  inserted  through  the  prcpellant  sample  as  illustrated  in  Figure  7. 
This  mounting  techrique  was  of  considerable  interest,  since  it  was  the 
favored  approach  in  the  ongoing  contractual  effort  at  the  time  this  analysis 
was  undertaken. 

The  basic  approach  followed  in  this  subsection  is  to  consider  the  pins  as 
beams  embedded  in  a semi-infinite  elastic  foundation  (the  propellant)  and 
then  to  examine  whether  the  deflections  thus  calculated  are  additive  to  the 
deflections  of  the  sample  considered  as  s beam. 

3.2.1  Stiffness  of  Foundation. 

Consider  the  mounting  pins  to  act  as  beams  supported  by  (embedded  in)  an 
elastic  foundation  of  infinite  extent,  as  shown  in  Figure  8(b).  A recessary 
parameter  for  analyzing  the  deflections  of  such  a beam  is  the  foundation 
modulus  k = q/w,  where  q is  the  load  per  linear  inch  exerted  on  the  found- 
ation by  the  beam  and  w is  the  foundation  deformation  produced  by  this  load. 
For  the  case  of  a uniform  bearing  pressure  p exerted  over  a square  area 
having  sides  of  length  2b  at  the  surface  of  a semi  - inf  in i te  uniform  clastic 
body,  the  average  deflection  over  the  loaded  area  is  (Reference  2,  p.  323) 

- = 1.90  pb  (l-v2). 

semi-inf  E 

For  the  case  cf  a pressure  p exerted  on  a plane  surface  at  the  center  cf 
an  infinite  uniform  elastic  body,  an  attractive  assumption  is  that  the 
pressure  load  would  be  equally  divided  between  the  two  semi  - inf  in i te  bodies 


comprising  the  Infinite  body  and  bounded  by  a plane  through  the  loaded  area. 
With  this  assumption,  the  net  result  for  the  infinite  body  is 


“l„f  ' [~ 


90  Pb  (l-u1 2)l  _ 0.95  pb  (l-v2) 


w 0.95  (IV) 

Finally,  recognizing  that  pb  - q and  taking  v - 1/2  and  E « Ef  where  Ef  is 
foundation  elastic  modulus. 


0.95  [l-(l/2n 


l.i*0A  Ef . 


2 

Note  that  k has  the  same  units  as  Ef  (lb/in  ) and  is  independent  of  the  size 
of  the  loaded  area  (i.e.,  of  b) . Actually  the  preceding  approach  to  deter- 
mining the  value  of  k was  replaced  by  another  approach,  discussed  in  3. 2. k. 
The  value  of  k actually  used  in  the  analyses  was 

k = 1.1»94  Ef . (26) 

3-2.2  Deflection  of  Pins  Inside  the  Propellant  Sample. 

The  pin  is  modeled  by  dividing  it  into  two  separate  portions,  the  "imbedded" 
portion  (Figure  8(b))  and  the  "non- imbedded"  portion  (Figure  8(c)).  The 
"imbedded"  portion  is  analyzed  as  a beam  on  an  elastic  foundation.  The 
assumptions  in  this  treatment,  as  discussed  by  Den  Hartog  (Reference  5, 
pp.  I4*»-11*5)  are  as  follows: 


1.  As  discussed  in  Reference  3.  page  58,  the  shearing  tractions  in  the  plane 
of  the  loaded  surface  vanish  for  the  case  of  v = 1/2.  Vanishing  of  the  shear 

tractions  in  this  plane  makes  the  indicated  assumption  valid  (see  Reference  4, 
pages  392  and  399) • 


(1)  Deflections  are  small. 

(2)  Plane  sections  remain  plane  under  deflection  (for  the  beam). 

(3)  The  foundation's  resisting  force  at  a point  on  the  beam  depends 
only  upon  the  transverse  deflection  of  the  beam  at  that  point. 

Den  Hartog  points  out  that  assumption  (3)  is  violated  by  the  behavior  of  a 
continuous  elastic  foundation,  but  states  that  the  excellent  agreement  of 
the  theory  with  experience  shows  that  the  assumption  produces  no  significant 


For  a beam  on  an  elastic  foundation,  it  is  convenient  to  work  with  the  non- 
dimensional  distance  Bx  where  \ is  the  physical  distance  measured  along  the 
beam's  length  and  B is  defined  by  (Reference  5>  page  143) 

r <”> 

where  k is  the  "foundation  modulus,"  E is  the  beam  material's  elastic  modulus, 
and  I is  the  moment  of  inertia  of  the  beam's  cross  section.  The  basic  theory 
treats  of  semi- inf  in i te  and  infinite  length.  If  the  beam  has  a finite  length 
h,  then  if  Bb  is  sufficiently  large,  the  semi  - inf  in i te  theory  holds  for  small 
values  of  6x1  that  is,  for  the  sem  i - inf i n i te  results  to  accurately  model  the 
deflection  and  slope  of  the  left  end  of  the  beam  in  Figure  8(b),  the  value  of 
Bh  must  be  "large."  Conversely,  if  3t  is  "small,"  the  beam  will  act  essen- 
tially as  a rigid  body,  with  its  displacement  governed  by  the  resisting 
forces  in  the  foundation.  Den  Hartog  provides  the  following  criteria  for 
"large"  and  "small"  values  of  6t  (Reference  5,  p.  161). 

Bl  large:  Bh  > 3 (28) 

Bt  smal I : Bh  < I . (29) 

The  two  types  of  behavior  are  discussed  separately  below. 
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3.2.2. 1 De t loot i on  of  Semi- I n f i n i te  Beam . 


i 


■ 


Since  the  load  P is  actually  applied  to  the  portion  of  the  pin  that  extends 
beyond  the  propellant  surface  (Figure  8(a)),  the  imbedded  portion  of  the 
beam  is  loaded  by  the  force  P and  the  moment  M - Pe,  as  shown  in  Figure  8(b). 
Superposing  the  results  given  by  Den  Hartog  for  P and  M treated  as  individual 
loads,  the  deflection  and  slope  at  any  point  x are 

[s,  ■ '<,<»*>  'T*  F3  (8*>]|x  (30) 

[V  T1  rl<s»)  ‘TT  "F1  <6*>]|x  01) 

with  the  values  ot  F j , F^,  and  F^  given  as  tabular  functions  of  8x  (Reference 
5,  page  146).  For  x = 0.  8x  " 0,  and,  as  it  turns  out,  F( (0x)  “ F2 (8x)  = 

F^ (8x)  ” I at  6x  = 0. 

Let  us  examine  the  effect  of  the  gage  standoff  distance  e on  the  pin's  deflected 

shape.  Consider  the  normalized  deflection  6 /#..  Making  the  proper  substitu- 

X “ 

tions  into  Equation  (30), 

Fi,(8x)  + jr-  f3(0x) 


or,  since  M - Pe, 

cSx  F/,(ex)  + F3(0x) 

«A  = I + Be 

Recall  that  the  beam  has  thus  far  been  regarded  as  semi-infinite.  It  is 
enlightening  to  choose  the  particular  value  x “ h and  examine  the  effects 
of  variations  in  8 and  e on  (6^/6^)  | This  has  been  done  in  Figure  9. 

utilizing  Den  Hartog's  tabulated  F.  values.  Before  discussing  Figure  9, 
however,  some  numerical  computations  are  in  order.  The  pins  used  at  Thiokol 
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in  the  Improved  Properties  program  were  #17  steel  pins,  having  a diameter 
d » 0.037  inch  (Reference  6).  Using  a modulus  of  elasticity  of  28  x 10^  psi, 

El  - E (55-)  - 2.576  lb-in2. 


Substitution  of  this  result  and  Equation  (26)  into  Equation  (27)  yields 


V'-^E  k 

6-v  trrw"  °-617  H 


Fairly  reasonable  values  for  maximum  propellant  modulus  E^.  and  gage  standoff 
distance  e are  10,000  psi  and  0.050  inches,  respectively.  These  values 
together  with  Equation  (33)  result  in  a value  of  Be  = 0.309.  For  the  same 
value  of  e and  a minimum  limit  of  the  modulus  of  100  psi,  Be  = 0.098.  For 
convenience,  the  values  Be  = 0 and  Be  = 0.5  were  considered  as  absolute 
upper  and  lower  bounds.  The  maximum  propellant  thickness  considered  for 
testing  was  0.5  inches.  Using  this  value  for  h and  making  use  of  Equation  (33), 


and 


Ef  Bh=l 

110  psi 

(34(a)) 

Ef  I Bh=3 

8,938  psi 

(34(b)) 

which  indicates  that  the  lower  bound  of  the  propellant  modulus  (100  psi) 
corresponds  to  a "short"  beam  and  rigid  beam  behavior,  while  the  upper  bound 
(10,000  psi)  corresponds  to  a "long"  beam  and  flexible  (semi- inf inite)  beam 
behav ior. 


Sufficient  background  has  been  developed  for  a discussion  of  Figure  9.  It 
is  evident  on  examination  of  Figure  9 that  the  relative  shape  of  the  deflec- 
tion curve  of  a sem i - inf i n i te  beam  is  weakly  influenced  by  changes  in  Be 
within  the  expected  range.  Figure  9 also  tends  to  corroborate  the  "long" 
and  "short"  beam  criteria  for  Bh  (6h  > 3 and  Bh  < I , respectively).  In 

-3/,- 


particular  (following  Den  Hartog's  rationale)  for  Bx  y 3,  bending  deflections 
(hence  bending  moment  and  shear)  damp  out  rapidly,  and  if  the  beam  were 
truncated  at  a point  in  this  region,  the  bending  moments  and  shear  forces 
thus  "destroyed"  would  have  a negligible  influence  on  the  beam  deflection 
curve  for  small  values  of  Bx;  then  the  deflection  curve  of  Figure  9 would 
apply  with  good  accuracy  to  the  finite  beam.  On  the  other  hand,  it  is  obvious 
that  truncation  of  the  beam  at  Bx  < I would  drastically  change  the  constraints 
on  the  left-hand  portion  of  the  bean;  both  bending  moment  and  shear  force 
would  be  significant  in  this  region.  For  the  latter  case,  the  deflections 
of  the  finite  beam  are  dominated  by  the  foundation  stiffness  rather  than  the 
beam  stiffness,  and  the  beam  itself  acts  essentially  as  a rigid  body  (Refer- 
ence 5,  p-  161). 

The  foregoing  discussion  indicates  that  there  are  two  rather  clear-cut  regions 
in  Bx  for  which  distinctly  different  models  apply.  Unfortunately,  the  area  of 
interest  is  the  muddy  region  in  between.  Den  Hartog  points  out  (Reference  5> 
p.  161)  that  for  this  region,  "In  many  practical  cases  a good  approximation 
is  obtained  by  using  either  the  stiff-beam  theory  (Bh  < 1)  or  the  long-beam 
theory  (Bh  > 3)  or  by  interpolating  in  between."  Further  discussion  of  this 
point  is  deferred  to  a later  part  of  this  subsection  of  the  report. 


3 . 2 . 2 . 2 Deflection  of  Short  Beam. 


Since  a "short"  beam  (Bh  < 1)  is  considered  to  act  as  a rigid  body,  its 
analysis  is  simple.  Given  that  the  foundation  obeys  Equation  (26),  the 
reaction  load  of  the  foundation  due  to  a beam  deflection  at  a point  is  pro- 
portional to  the  deflection.  The  analysis  is  most  easily  done  by  applying 
the  superposition  principle  as  in  Figure  10(a).  Static  equilibrium  equations 
for  the  individual  free-body  diagrams  of  Figure  10(b)  and  10(c)  result  in 
the  following  equations: 


9,  h 


k6. 


EFy  = P 

"o  - 2 f“\  (<  - ?)“  (*  -t)-  »c 


EM 


k«2  h 
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Figure  10.  Analysis  of  Rigid  Beam  on  Elastic  Foundation: 

(a)  Decomposition  of  Load  into  Central  Force 
and  Moment; 

(b)  Free  Body  for  Beam  Loaded  by  Force  P; 

(c)  Free  Body  for  Beam  Loaded  by  Moment  M . 


by  combining  the  deflections  of  the  imbedded  and  nonimbedded  portions  of  the 
beam.  The  imbedded  portion  has  been  treated  previously;  the  nonimbedded 
portion  will  now  be  discussed. 

3.2.3* I Deflection  of  Nonimbedded  Portion  of  the  Pin. 

Referring  to  Figure  8(c),  let  6'^,  be  the  deflection  at  point  A’  for  a beam 
cantilevered  from  point  A.  For  the  loading  shown, 

Pe^ 

5 A'  = 2EI  (37) 

3-2. 3.2  Combined  Results  for  the  Entire  Pin. 

In  3.2.3. 1 above,  the  nonimbedded  portion  of  the  pin  was  considered  as  a 
cantilever  beam.  Since  the  beam  is  actually  connected  to  the  imbedded  por- 
tion of  the  pin,  the  vertical  motion  («A)  and  rotation  (eA)  of  point  A must 
be  incorporated.  The  final  result  for  the  deflection  at  point  A'  is 

Pe2 

6A'  = 2ET+  6A  e + 6A  (38) 

For  the  short  (rigid)  beam  case,  6^  and  0^  are  calculated  from  Equations  (35) 
and  (36).  For  the  long  beam  case,  6 A and  6 A are  calculated  from  Equations  (30) 
and  (31)  with  x - 0 (as  noted  earlier,  this  results  in  Fj  (ex)  - F2  (bx)  = 

F3  (3x)  = fj,  (BX)  - I).  The  complete  results  are  summarized  in  Equations  (39) 
and  (AO) : 


Long  Beam  Case  ("infinite  beam"): 


6A' 


Pe2  . ( 202h2P  . Af33h2Pe 
inf  2EI 


kh 


) e * + 26h2Pe\ 

I \ kh  kh2  I 


(39) 


Short  Beam  Case  (rigid  beam): 


V 


rigid 


Pe2  P_ 
2EI  kh 


P_  /6  12e\  ^ P /.  6e\ 

kh  lh  77  ^ ( h“) 


(AO) 
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For  intermediate  cases,  i.e.,  cases  for  which  I <_  f?h  <_  3,  the  approach  (as 
mentioned  earlier)  suggested  by  Den  Hartog  is  to  interpolate  between  these 
two  results.  Applying  this  advice. 


Also, 


and 


where 


V “ Kdh  (6A'  , ' V . .}  + V . 

inf  rigid  rigid 


«A  " *Bh  UA  ■ SA  ) + SA 
* dh  Alnf  Arigid  Ar  ig  id 

% “ Kdh  <Y  f - eA  . . > + °A  . .. 

inf  rigid  rigid 


lSB  " Kdh  ^B.  , ‘ lSB  . .J  + 6B  . . . 

inf  rigid  rigid 


(41) 

(42) 

(43) 

(44) 


0, 

(3» 


Kgh  =*  (ph-1 ) /2 , 


0 < 6h  < 1 

1 < Bh  < 3 
3 < 6h 


(45) 


Since  the  end  use  of  the  calculations  made  according  to  the  above  relation- 
ships is  to  determine  the  net  effect  on  the  total  strain  measurement  bias,  a 
useful  parameter  is  the  additive  component  of  strain  bias  in  inches/inch  per 
pound  of  force  exerted  on  the  pin  by  the  clip  strain  gage.  Since  there  are 
two  pins  which  deflect  (at  the  loading  point)  through  a distance  5^,,  the 
additive  bias  is  given  by 


- "DELTA  STRAIN  BIAS  t P" 
2 ( <S  A , ) / P 


(46) 


Where  Lq  is  the  total  distance  between  pins,  i.e.,  the  gage  length  of  the 
clip  gage.  (All  calculations  performed  used  a fixed  value  of  0.400  for  Lq.) 
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The  manipulations  required  to  combine  Equation  (46)  with  the  results  of  the 
earlier  analyses  of  sample  bending  deflection  to  calculate  the  total  strain 
measurement  bias  is  deferred  until  a later  subsection  of  this  report. 


Appropriate  calculations  to  evaluate  Equations  (41),  (42),  (43),  and  (46) 
were  programmed  on  a programmable  calculator  In  accordance  with  the  rela- 
tionships developed  above.  Selected  results  (e  - 0.050")  are  shown  in 
Figure  11.  It  can  be  seen  that  the  interpolated  result  for  transitions 
very  smoothly  from  the  rigid  beam  result  to  the  infinite  beam  result  as  Bh 
goes  from  I to  3.  The  behavior  of  this  model  of  the  steel  pin  in  the  pro- 
pellant specimen  appears,  in  a general  sense,  to  be  quite  satisfactory. 

Figure  12  provides  a physical  illustration  of  the  pin  deflection  as  described 
by  the  model.  This  figure  shows  the  predicted  vertical  deformation  of  the 
pin  treated  as  a "rigid  beam,"  an  "Infinite  beam,"  and  as  an  intermediate 
case  ("interpolated  result").  As  noted  on  the  figure,  the  particular  results 
shown  are  for  Ef  - 1000  psi;  the  calculated  values  of  Bh  and  Be  are,  respec- 
tively, 1.735  and  0.174.  (Actually,  the  calculator  program  provides  only 
the  deflections  at  points  A and  A'  and  the  slopes  at  point  A for  the  three 
deflection  curves.  The  remainder  of  Figure  12  was  constructed  by  interpolating 
6/6^  from  Figure  9 for  Be  » 0.174  to  obtain  the  infinite-body  curve,  then 
manually  interpolating  between  this  curve  and  the  rigid-body  "curve"  for 
Bh  = 1.735.)  The  basis  of  Den  Hartog's  discussion  of  the  effects  of  removing 
the  part  of  the  beam  which  lies  outside  the  foundation  (see  3-2.2. 1)  is 
clearly  evident  in  Figure  12. 


3.2.4  Re-examination  of  Foundation  Stiffness. 

The  earlier  discussion  of  foundation  stiffness  (3-2.2)  cited  a very  simple 
expression  for  the  foundation  stiffness.  This  subsection  concerns  itself 
with  severa1  factors  which  were  examined  in  detail  to  obtain  a final  estimate 
of  the  effective  foundation  stiffness  (a  particularly  important  factor  being 
the  interference  between  two  opposed  points  or  distributed  loads  acting  within 
the  medium  or  foundation).  This  subsection  also  examines  the  manner  in  which 
the  sample  bending  deformation  and  the  localized  deformations  of  the  pin  and 
surrounding  propellant  interact  to  produce  the  total  motion  of  the  pin. 
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3. 2. *4.1  Deflection  Under  Highly  Localized  Loads. 


The  deflection,  produced  by  a concentrated  (point)  load  P acting  at  the 
surface  of  a semi- inf  ini te  medium,  which  exists  at  a particular  point  in 
the  medium  (r,  z,  9)  is  given  by  (Reference  4,  p.  402): 


“ ' 2"£f  [ 


? 2 2 2 2 2 
(1  + v ) z (r  + z)  + 2 (1  - v ) ( r + zL)  (47) 


where  the  force  acts  in  the  positive  z direction  at  the  origin  of  coordinates 
and  the  medium  has  an  elastic  modulus  E and  a Poisson's  ratio  v.  Consider 
the  deflection  (see  Figure  13(a))  of  a point  on  a hemispherical  surface  of 
radius  a,  centered  at  r = z = 0.  For  Poisson's  ratio  = 0.5,  it  can  be  shown 
that  Equation  (47)  results  in 


“ surf  ' fel7  [?  <‘»s2  » * "] 


Two  quantities  are  of  interest:  w jsurf  or  the  deflection  averaged  over  the 
hemispherical  surface;  and  w arc  or  the  deflection  averaged  over  an  arc  in 
a plane  perpendicular  to  the  hemisphere's  base.  Integration  of  Equation  (48) 
(with  the  proper  substitutions)  results  in 


rr\ 


ad<f>d© 


whence 


rv/1 

1 w 


" I surf 
na/2 


1.125  w 


The  average  deflection  for  a circular  area  of  radius  a loaded  by  a uniform 
pressure  p is  given  as  (Reference  4,  page  409): 

0.54  P ( 1 — y2) 
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(a)  COMPUTE  AVERAGE  DEFLECTION  OVER  A CIRCULAR  ARC 

(RADIUS  a=  PIN  OIA  2)  PRODUCED  BY  A CONCENTRATED 
FORCE  P. 


[ 


J 


(b)  COMPUTE  AVERAGE  DEFLECTION  OVER  CIRCULAR  AREA  OF 
RADIUS  a UNDER  UNIFORM  PRESSURE  p=P/(*a2' 


I* 


(c)  CALCULATE  DEFLECTION  UNDER  A RIGID  DIE  OF  RADIUS 
a LOADED  BY  A CONCENTRIC  FORCE  P. 


I 


Figure  13.  Approaches  to  Determination  of 
Foundation  Stiffness  Under  a 
Highly  Localized  Load 


This  result  applies  to  Figure  13(b). 


Use  of  Equations  (49)  and  (50)  result  in: 


,/w  = 0.786 

surf/  press 


and 


/w  = 0.884. 
arc/  press 


The  deflection  for  the  rigid  die  (Figure  13(c)) 
(Reference  2,  page  323): 


= P P-v2) 
Vigid  2aE, 


(51) 

(52) 

under  a load  P is  given  by 


Again  using  Equations  (49)  and  (50)  and  taking  v = 1/2, 


surf/  rigid 


0.849 


0.955 


(53) 

(54) 


It  appears  that,  for  the  same  total  force,  the  technique  of  averaging  the 
deflections  due  to  a point  load  over  the  arc  of  radius  a produces  nearly  the 
same  result  as  obtained  for  the  deflection  of  a rigid  die  of  the  same  radius 
and  a result  close  to  that  for  a circular  area  of  radius  a under  a uniform 
pressure. 


Timoshenko  and  Goodier  show  that  Equation  (47),  when  integrated  for  the  case 
of  Figure  13(b),  produces  the  following  results  along  the  loaded  surface 
(Reference  4,  pages  403  - 409) 


(1-v2) 


jja_  2 (1-v  ) P 


/ r*  1-  \ 


iiifc. . — dBH 


Mi- 


[F'  & r)  - ('  • 7)  (p  f)] 


for  a < r < °°,  and 


MM  [F'  ®] 


for  0 < r < a 


where  F|  denotes  an  elliptic  integral  of  the  fi-st  kind  and  F^  denotes  an 
elliptic  integral  of  the  second  kind.  Equations  57  and  58  were  normalized 
wi  th respect  to  w (Equat ion (55))  and  the  values  of  this  normalized  deflec- 
tion were  plotted  to  produce  Figure  14.  The  deflection  produced  by  an 
equivalent  point  load  (Equation  (47))  is  also  shown  in  Figure  14.  It  is 
noticeable  that  the  two  curves  are  almost  indistinguishable  for  r/a  greater 
than  3.  Another  noticeable  feature  is  the  rapidity  with  which  the  local 
deflections  die  out.  Before  drawing  any  conclusions,  however,  the  effects 
of  an  opposing  load,  coaxial  with  the  load  of  interest  but  acting  at  a 
distance  Lq  along  the  axis,  must  be  considered.  As  a first  step,  the  deflec- 
tion at  points  on  the  z axis  (see  inset.  Figure  14)  can  be  found  by  integrat- 
ing Equation  (47),  with  the  result 


I7  (i 


/ 2 ^ 2 
a + z 


O-v2)^ 


/ 2 ^ 2 

a + z - z 


It  is  easily  demonstrated  that  this  equation  produces  the  same  value  of 
w = w for  r = 0,  z = 0 as  is  given  by  Equation  (55).  Taking  v = 1/2  and 
normalizing  with  respect  to  w , 


(w/"«ax)r,0  ' V ’ * <t) 


V'  * (f) 


Equation  (60)  is  plotted  in  Figure  15,  along  with  the  normalized  deflection 
produced  along  the  z axis  by  an  equivalent  point  load.  Note  that,  as  in 
Figure  14,  the  point  and  distributed  loads  produce  essentially  identical 
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of  Symmetry  Produced 


results  at  a characteristic  dimension  “ = 3 and  beyond.  Figure  15  also 
shows  the  results  of  superposing  the  deflections  for  two  distributed  loads: 
one  acting  downward  at  z - 0 and  one  acting  upward  at  z = Lq  = 0.4,  or 
z/a  = LQ/a  - 0. VO. 0185  - 21.62,  i.e.: 


net 


r-0 


r-0 

z-0 


w 

w 

max 


z-0 


Comparing  the  deflection  curve  for  the  two  opposed  load  distributions  with 
that  for  the  single  load  distribution,  it  appears  that  the  deflection  under 
each  load  is  decreased  by  the  presence  of  the  opposing  load.  The  effective 
stiffness  of  the  medium,  in  other  words,  is  increased  by  the  interference 
between  the  two  loads. 

Taking  advantage  of  the  rapid  convergence  between  the  distributed-load  and 
point-load  deflection  curves  as  z/a  increases  (Figure  15).  the  distributed- 
load  deflection  at  the  loaded  surface  (Figure  14)  was  corrected  for  the 
presence  of  an  opposing  distributed  load  at  z/a  = LQ/a  by  substituting  an 
equivalent  point  load  for  the  distributed  opposing  load.  The  result  is 
shown  in  Figure  16.  Also  shown  in  Figure  16  is  the  normalized  average 
deflection  of  the  uniaxial  specimen,  computed  by  integrating  the  strain 

given  by  Equation  (2)  over  the  length  Lq.  (P1  = 2P  was  substituted  for  P 

in  Equation  (2)  to  account  for  the  use  of  a semi-infinite  body  model  rather 
than  an  infinite  body  model  in  the  above  analyses.  Notice  also  that 
Figure  16  attributes  the  final  result  to  the  infinite  body  rather  than  to 
the  semi -inf ini te  body.  As  discussed  in  3.2.1,  this  procedure  is  valid 
only  for  a Poisson's  ratio  of  1/2.) 

Comparing  Figures  15  and  16,  it  can  be  seen  that  the  average  deflection  over 

the  loaded  area  is  reduced  from  0.85  w to  0.799  w by  the  interference 

o o 

effect  of  the  second  load.  This  represents  an  increase  of  6 percent  in  the 
effective  stiffness  of  the  medium.  It  is  also  apparent  that  the  interference 
from  the  second  load  causes  the  deflection  curve  to  "damp  out"  somewhat  more 


rapidly  as  r/a  increases.  These  results  indicate  that  the  effects  of  this 
interference  must  be  considered  In  the  analyses  of  the  mounting  pins. 


Figure  16  provides  a basis  for  assessing  the  importance  of  "edge  effects" 
on  the  validity  of  the  infinite  medium  model  for  analyzing  the  interaction 
of  the  pin  with  the  finite  sample.  It  is  evident  in  Figure  16  that  the 
deflection  is  negligible  for  values  of  r/a  greater  than  5.  We  can  therefore 
infer,  in  an  argument  parallel  to  the  one  Den  Hartog  uses  for  the  infinite 
beam  (see  3.2.2),  that  material  beyond  r/a  = 5 could  be  physically  removed 
with  little  effect  on  the  average  deflection.  In  fact,  it  would  appear 
qualitatively  that  removal  of  all  material  to  the  right  of  a plane  tangent 
to  r/a  of,  say,  2 or  3 would  still  have  little  effect  on  the  average  deflec- 
tion or  effective  medium  stiffness.  However,  the  impact  of  removing  material 
would  rapidly  increase  for  smaller  values  of  r/a  and  would  certainly  be 
significant  as  the  plane  approached  r/a  = 1.  We  can  therefore  conclude  that 
the  infinite-medium  result  applies  accurately  for  a localized  load  provided 
the  loaded  area  is  at  least  two  to  three  times  its  own  radius  distant  from 
a free  surface,  but  that  the  infinite-medium  result  underpredicts  the  deflec- 
tions when  the  loaded  area  is  closer  to  the  surface.  This  conclusion  will 
be  of  use  in  the  following  subsection  of  this  report. 

One  final  useful  observation  can  be  made  based  on  Figure  16.  The  normalized 
deflection  curve  depicts  the  local  deflections  ofan  infinite  medium  under  a 
highly  localized  load.  However,  nearly  identical  localized  deflections  must 
occur  within  a finite  sample  under  the  same  load,  provided  it  is  sufficiently 
large  that  edge  effects  can  be  ignored  (as  discussed  above).  These  conditions 
are  satisfied  for  the  1/2-inch  by  1/2-inch  propellant  sample,  as  is  evident 
from  Figure  16.  Based  on  St.  Venant's  principle  (Reference  5,  page  117)  we 
would  expect  the  local ized  effects  of  the  concentrated  load  to  be  merely  a 
perturbation  on  the  average  sample  behavior.  In  other  words,  if  we  consider 
the  deflection  in  the  finite  sample  to  be  given  by  the  sum  of  two  components. 


w wuniaxial  + wlocal 


(61) 


J 
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where  w , . . is  given  by  Equation  (2),  we  would  expect  that 

Un  I dX I a I 


w = w . . . = w uniaxial  + w.  . 

uniaxial  local 


or 


w,  , = 0. 
local 


However,  numerically  averaging  the  deflection  curve  of  Figure  16  produces 
the  result 


w/w 


0.0349 


r=0. 25 

while  the  uniaxial  deflection  for  the  1/2  x 1/2  sample  was  computed  as 


w . . . /w 

uniaxial/  max 


0.056 


z=0 


It  appears,  then,  that  the  deflection  calculated  for  the  infinite  medium  is 
more  nearly  equal  co  the  total  local  deflection  of  the  sample  than  to  the 
component  w|oca|‘  ln  fact,  the  component  of  deflection  w|oca|  's  approxi- 


loca 

mate  I y given  by 


w = w ^ w 

local  inf  medium  uniaxial 


This  result  will  be  utilized  in  the  following  subsection. 


(62) 


i 


3 . 2 . 4 . 2 Deflection  Under  Distributed  Localized  Loads. 

Since  the  major  part  of  both  the  localized  and  the  average  deformation  of 
the  sample  is  in  flexure  rather  than  axial  deformation,  a foundation  stiff- 
ness derived  from  response  to  a distributed  bending  moment  appears  to  be 
appropriate.  The  simplest  load  of  this  type  is  the  triangularly-distributed 
line  load  of  Figure  17(a).  Rewriting  Equation  (47)  (in  the  differential  form) 
in  terms  of  the  variables  shown  In  Figure  17  and  substitution  of  v ■=  1/2 
results  in  the  differential  equation  for  deflection  within  a semi  - inf  in ite 
med i urn: 


V 


1 
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(a)  Triangular  Load  Distribution  and  Geometry 


y-y/h 


ds-adP 


J 

Z=Z/h 

(b)  Integration  Contour  for  Calculation  of  w„ 


Figure  17.  Nomenclature  for  the  Analysis  of  Deflec- 
tions Under  a Triangular  Load  Distribution 


r 


1.5  q, 


-3/2 


-1/2- 


dw 


nEh 


max  'z2  (A2  + y2  + z2)  + (A2  + y2  + z2)  (A+x)  dA. 


\f 2 2 

Further  substitution  of  a *vy  + z , x " x/h»  Y “ y/h,  z “ z/h,  a ■*  a/h  and 
A ■ A/h  produces  the  equation 

1.5  q__v  r*,  "3/2  -1/2-1  « 

dw  “ ^ z (A  + a ) + (A  + a ) (A  + x)  dA  (63) 

Integration  of  Equation  (63)  over  the  length  of  the  clip  gage  mounting  pin, 
i.e.,  from  A(  «*  -(x  + j)  to  A2  * (y  - x).  produces  the  following  result: 


1-5  q, 


rrE 


max  ) '2 
< z 


• 

H - 

+ Z2  1L 

- 

^2 

+ 

" 

-1 

A 

VA2  ♦ a2 

V':2  + "2 

A a 

2 

a 

V'2  ^ "2 

A + a 

L J 

L 

A, 

— — 

X *n 


If  we  now  define  w = w 
o 


in  (a  +VA2  ♦ 


(64) 


Jxl 


X - 1/2 

a - a . = a 

pin/h  o 


i.e.,  the  value  of 


w at  the  end  of  the  pin  and  lying  on  a surface  in  the  propellant  correspond- 
ing to  the  pin's  surface,  it  can  be  shown  from  Equation  (64)  that 


1.5  q, 


irE 


ma*  [(aQ2  sin2e)  ki  + k2 ] = k [(aQ2  sin2°^  ki  + ka]  ^5) 


where  0 is  defined  in  Figure  17(b).  Integrating  Equation  (65)  with  respect 

to  ds  ■ aQd0  from  8 = 0 to  0 ■ n/2  and  dividing  by  (ira/2)  yields  the  following 

average  value  for  w : 

o 


w * k (2a  ) 
o 
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i j 
1 

1 

A 


Evaluation  of  kj  and  k£  from  Equation  (6A)  with  a = (.01 85) / (0 . 5)  ■ 0.037, 
Aj  » -I,  and  A,  » 0 yields  the  result 

w = 1.2635  k - 1.2635 
o 

and  Equation  (61*)  becomes 


w/wq  - 0.7915  < z 


- 

x 

- 

-1 

i ...  L _ 

V'2  '2 

a2  Vr2  *2 

A + a 

x. 

A + a 

T X„ 


♦ K2  + ;2 


r v^2  *2\1  2 I 

x In  [a  +vA  + a ) A > 
- ^ 1 ' 


(66) 


Equation  (66)  was  programmed  on  a programmable  calculator  and  calculations 
were  performed  both  for  a single  triangular  load  distribution  and  for  two 
opposed  triangular  load  distributions  separated  along  the  z-axis  by  a distance 
of  * 0.4  inches.  The  results  are  shown  in  Figures  1 8 , 19,  and  20.  The 
influence  of  the  second  line  load  (compare  Figure  18(a)  with  Figure  18(b))  is 
quite  similar  to  that  found  for  interacting  highly  localized  loads  in  the 
previous  subsection;  i.e.,  deflections  are  slightly  reduced.  Figure  19  shows 
the  variation  in  deflection  along  the  (imaginary)  "pin  surface"  at  the 
extreme  end  of  the  pin  axis.  There  is,  of  course,  significant  distortion  of 
this  surface  as  the  medium  deflects  under  the  line  load,  although  this  dis- 
tortion is  considerably  less  than  the  distortion  of  a circular  arc  around  a 
point  load  (this  comparison  is  shown  in  Figure  19).  Since  it  was  shown  in 
the  preceding  subsection  that  the  average  deflection  over  the  circular  arc 
produced  by  the  point  load  closely  duplicated  the  deflection  under  an  equiva- 
lent ly-loaded  rigid  die,  averaging  deflections  over  this  imaginary  "pin  surface" 
should  provide  a good  estimate  of  the  deflection  under  a rigid  die  of  width 
2a  and  length  h,  subjected  to  a load  equivalent  to  that  in  Figure  17(a).  The 
calculated  deflections  along  the  "pin  surface"  for  four  values  of  6 from 
0°  to  90°  are  plotted  in  Figure  20.  As  noted  from  Figure  19,  there  is  distor- 


DEFLECTION  for 
SAMPLE  FLEXURE 


e Triangul 


ons  at  z»l 


DEFLECTION  FOR 


SAMPLE  FLEXURE. 


(b)  Deflect ions 


Opposed  Triangular  Distributions 


Figure  18 


Deflect  ions  Produced  in 
the  .Triangular  Load  Dt s' 


Fi  gure  I 7^ ( a 


Hi] 

jljj 

w/w  FOR  POINT  LOAD 
(FROM  Eq ( 48)  ) 


TWO  OPPOSING 
LINE  LOADS 


0,  DEGREES 


Dtstri&MtfoH-Atewtg  an  Arc  of  Rod 


>i)i  Produced 


He  Triangular  Load  Divtr 


lion  of  figure' 


EXTRAPOLATED  ROTATION 


AVERAGE  ROTATION 


DEFLECTION  FOR 
SAMPLE  FLEXURE 


0 

i i 

0.1  0.2  0. 

X-X/h 

1 

3 0.4 

i 

0.5 

Figure  20 

. Deflections  on  a Cyli 
0.037  Produced  by 
bution  of  Figure  1 7(a 

ndrical  Surface  o 
the  Trianqular  Lo 
). 

f 

ac 

Rad i us 
Distri 

; 

j . : ..  ■ 

| 

. j 

j.  • i ; 

; . . j ... 

. ' :• 

for  the  Intended 


Th I a La  1 


tion  of  each  cross  section  of  the  "pin"  as  it  deflects.  There  is  also  a 
substantial  bending  deformation  which  would  not  occur  with  a rigid  die. 

One  approach  to  the  desired  stiffness  calculation  is  to  calculate  the  average 
rotation  over  the  length  of  the  "pin"  and  take  this  value  as  the  "rigid  die" 
deformation.  This  calculation  was  performed  by  averaging  the  four  deformations 
at  each  cross  section,  computing  (w/wq)/x,  and  averaging  these  values  to 
compute  a rotation.  In  the  course  of  these  calculations,  it  was  found  that 
the  deflected  shapes  of  the  cross  sections  were  very  nearly  geometrically 
similar:  (w/wq) /\  had  almost  exactly  the  same  shape  at  each  cross  section. 

While  the  above  approach  (which  yields  the  line  labeled  "average  rotation" 
in  Figure  20)  is  straightforward,  it  neglects  the  "edge  effects"  discussed 
in  the  preceding  subsection  of  this  report.  As  noted  in  the  discussion  of 
Figure  16,  these  effects  begin  to  be  important  as  the  distance  from  the  free 
surface  decreases  below  two  to  three  times  the  radius  of  the  loaded  area. 

The  corresponding  values  of  x in  Figure  20  are 


- - 2 
a 


0.5  - 2 (.037)  = 0.1*26 


0.5  - 3 (.037)  = 0.389 


It  appears,  therefore,  that  the  "edge  effects"  would  begin  to  reduce  the 
localized  stiffness  of  the  medium  (i.e.,  to  increase  the  localized  deflections) 
at  approximately  the  point  for  which  the  deflection  curves  of  Figure  20  begin 
to  "bend  over."  Furthermore,  the  localized  effective  stiffness  would  be 
expected  to  drop  to  about  half  the  value  in  the  center  of  the  sample  at  the 
edge  of  the  sample  (i.e.,  the  deflections  would  double).  The  net  effect  would 
be  to  approximately  negate  the  departure  of  the  deflection  curves  in  Figure  20 
for  x/h  greater  than  3.  Based  on  this  reasoning,  and  on  the  observed  similarity 
in  the  cross-sectional  shapes,  the  maximum  deflection  for  stiffness  calculation 
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purposes  was  determined  as  follows:' 


(1)  From  Figure  20,  It  was  determined  that  w * w when  6 = 65° 

o o 

(2)  This  value  of  9 was  used  to  calculate  (w/w  ) 

O X B U • z 

(3)  This  value  of  w/wq  was  then  used  to  extrapolate  w/wq  to  x ■ 0.5, 

obtaining  (w/w  ) ■ 2.22 

3 o max 

The  value  thus  obtained  Is  considerably  hl*ner  than  the  "average  rotation" 
value  of  1.90  (see  Figure  20). 

3 • 2 . 4 . 3 Effective  Foundation  Stiffness. 

The  preceding  subsection  provides  a calculated  deflection  distribution  pro- 
duced in  a semi-infinite  medium  by  the  triangular  load  distribution  of 
Figure  17(a).  If  we  take  the  deflection  as  an  effective  linear  distribution, 
the  foundation  stiffness  is  given  by  the  definition  of  k: 


k = q/w  = 


(w/w  ) (w  /k) k 
o o 


Where  the  factor  2 is  applied  to  q to  obtain  the  stiffness  for  an  infinite 

max 

medium  instead  of  a semi- inf inite  medium.  From  the  preceding  subsection, 
(w/wq)  = 2.20,  (wQ/k)  = 1.2635,  and  k = (1.5  %iax)/(1,^f)  • Substitution  of 
these  values  into  the  above  equation  produces  the  result 

k - I.W  Ec  (67) 


2.  In  retrospect,  step  (1)  of  this  procedure  was  based  on  an  erroneous 
assumption^  that  unsymmetrical  distortion  was  responsible  for  w not  being 
equal  to  wq  at  0 - 1*5°.  Actually,  as  is  evident  in  Figure  19.°the  shift  in 
w/wq  is  simply  the  result  of  Interference  from  the  opposing  line  load.  The 
resulting  value  of  (w/w_)max  's  a^out  10  percent  too  high  (a  small  error  in 
view  of  the  uncertainties  Involved). 
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For  comparison,  values  of  k can  be  calculated  for  the  highly  localized  loads 
analyzed  in  Define  the  load  per  unit  length  for  the  circular  areas 

of  Figure  13  as  q * P /l  where  P is  the  total  load  and  l Is  the  length  of  an 
equivalent  square  (i.e.,  a square  having  the  same  area  as  the  circle).  Then 


vMa) 


and,  for  the  average  deflection  produced  by  a point  load  acting  within  an 
infinite  medium,  we  have,  from  Equation  (50) 


w ” 2 semi  - inf i 


1 . 1 ( 9P_\  „ I (9 ^ n 

infinite  medium'  2 \8irE  ^a  / 2 \8^/ifE^  / q 


k * q/w 


2 (Vi)  Ef 


3-151  E. 


(point  load,  deflection  averaged  over  arc  of  radius  a). 

Similar  calculations  result  in  the  following  stiffness  values: 

Uniform  pressure  on  a circular  area:  k = 2.786  E^ 

Deflection  of  a rigid  circular  die:  k = 3.009  E^. 

The  selected  technique  (Equation  (67))  of  calculating  stiffness  of  the  medium 
therefore  results  in  a stiffness  about  half  that  for  a highly  localized  load 
on  a circular  area.  However,  the  selected  stiffness  is  very  close  to  the 
value  obtained  for  a square  area  loaded  by  a uniform  pressure. 

3 • 2 . . 4 Ad  justment  for  Sample  Flexural  Deformation. 

The  triangular  load  distribution  of  Figure  17(a)  is  statically  equivalent  to 
a pure  bending  moment  given  by 


* ' 


' 


fh/Z2  q q h2 

\’2J  ~1~  ^ 

actmg  on  the  semi  - inf  in  i te  body  of  the  analysis.  The  moment  acting  on  an 
infinite  body  (the  propellant  sample)  is  twice  this  value,  or 

q h2 

M - 2M  - -222 • 


From  the  analysis  of  3.1,  the  bending  moment  in  the  sample  gage  length  is, 
since  Rj  = 0,  (Equations  (3)  and  (5)) 


M = -M.  + R z + M = -M.  + M 
I I o I o 


where  (Equation  (4)) 
M 


"i  - - 7 h - '21 


substitution  of  the  above  expression  for  M yields 

o 

-2  M a -2  q h2a 
M o ^max 


31 


and  the  surface  bending  strain  in  the  center  section  of  the  sample  is  given 
by 


(Ae  ),  . 

z bending 


Mh 
2Ef  I 


'max 


h3a 


3Ef  It 


The  surface  axial  deflection  at  z = a (see  Figure  1)  is  given  by 

(w  ) = L°  (ae  ) = qmax  h a Lo 

max  sample  2 '“V  bending  STJjt  * 

where  a factor  of  -1  is  incorporated  because  a negative  strain  produces  a 
positive  value  of  w. 
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Recognizing  that  I » h / 1 2 , 


2 q a L 
max  o. 


max'  sample 


and,  since  wq  - 1.2635  (1.5  qmax)/(nEf), 


2 7i  a L 


L.  n a l. 

^ wmax ^ s«t.np  1 e ^Wo  J ( 1 . 2635)  ( l .5 


substituting  a - 1.80  inches,  Lq  - 0.4  inches,  h - 0.5  inches,  and  l • 4 inches 


produces  the  final  resu't 


[(w  ) . /w  1 - 1.193 

L max  sample  oj 


This  sample  bending  displacement  is  shown  in  Figures  18  and  20. 

If  we  apply  the  same  logic  as  in  3-2.4. 1 the  displacement  in  the  z direction 
is  given  by 


w = w . + w.  , 1°°' 

sample  local 

where  w , is  the  bending  deflection  of  the  sample  treated  as  a beam,  and 
sample 

wlocal  * winf  medium  wsample  ^winf  medium^  C 


where 


W — w 

inf  medium  sample 
winf  medium 


and,  substituting  computed  values  for  the  deflections, 
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3.2.5 


Added  Bias  in  Strain  Gage  Measurement. 


The  preceding  discussion  provides  the  groundwork  for  computing  the  effect  of 
local  deflections  on  the  measurement  bias  produced  by  clip  gage/sample  inter- 
action. The  bias  due  to  local  deflections  is  considered  as  an  additive  effect 
which  is  superimposed  on  the  bias  predicted  by  modeling  the  sample  as  a beam. 

As  shown  in  3. 2. 4. A,  the  local  component  of  deflection  of  the  sample  under  a 
triangular  line  load  is  considerably  less  than  the  deflection  of  an  infinite 
medium  under  the  same  load.  The  logical  application  of  this  result  would  be 
to  apply  the  calculated  reduction  factor  (Equation  (70))  to  the  pin  deflec- 
tions calculated  earlier  (in  3-2.3).  Examination  of  Figure  12,  however, 
reveals  a basic  conflict  with  this  concept.  The  rigid  beam  model,  as  shown 
in  Figure  12,  produces  a pin  deflection  curve  which  closely  resembles  the 
deflection  of  a beam  cross  section  under  bending  loads.  The  motion  of  the 
pin  is  simply  a magnified  version  of  the  cross-sectional  deformation.  However, 
it  is  obvious  in  Figure  12  that  this  observation  is  not  true  for  the  infinite 
beam  model:  the  beam  itself  undergoes  no  rigid  body  motion,  either  in  rotation 
or  vertical  deflection.  It  therefore  is  apparent  that  the  rigid  body  motion 
of  the  pin  produced  by  deformation  of  the  sample  as  a beam  must  be  added  to 
the  deflection  shown  for  the  infinite  beam  in  Figure  12  to  obtain  the  total 
motion  of  the  pin.  Conversely,  the  local  component  of  the  pin  deflection  is 
equal  to  the  deflection  given  by  Figure  12  without  the  reduction  factor  of 
Equation  (70). 

The  calculator  program  used  to  generate  Figure  11  was  modified  to  apply  the 
reduction  factor  of  Equation  (70)  to  the  "rigid  beam"  deflections  and  slopes; 
the  computed  deflections  and  slopes  for  the  infinite  beam  model  were  left 
unchanged.  The  deflections  and  slopes  were  then  interpolated  by  the  calculator 
program  in  the  same  way  as  in  Figure  11.  The  results  for  an  eccentricity  of 
0.050  inches  are  shown  in  Figures  21  and  22.  As  before,  the  interpolation 
yields  results  that  appear  reasonable. 

Figure  23  shows  the  bias  in  the  strain  measurement  produced  by  local  deforma- 
tions alone,  for  several  values  of  propellant  modulus,  as  a function  of  eccen- 
tricity. These  results  were  calculated  (by  the  calculator  •'-ogram  discussed 
above)  according  to  Equation  (46). 
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Figure  23.  Predicted  Strain  Bias  Caused  by  Local  Deflection 
of  Mounting  Pins  (h=*0.50  in,  L =0.40  in) 


Net  Gage-Sample  Interaction. 


In  3.1  and  3.2  of  this  report,  analyses  were  conducted  to  estimate  the 
separate  contributions  to  bias  in  clip  strain  gage  measurements  on  a uni- 
axial tensile  specimen  caused  by  sample  deformation  (as  a beam)  and  by 
local  deformation  of  the  clip  gage  mounting  pins.  The  analyses  in  the 
present  subsection  combine  these  separate  contributions  to  obtain  the 
net  bias  in  the  clip  gage  measurements.  The  appropriate  equations  are 
developed  and  calculations  are  carried  out  for  two  selected  gage  stiff- 
nesses. In  addition,  the  pin  deformation  and  displacement  resulting 
from  both  sample  deformation  as  a beam  and  localized  effects  is  compared 
with  three-dimensional  finite  element  analysis  results. 

3.3*1  Combined  Bias  In  Strain  Measurement. 


The  total  bias  in  the  strain  gage  measurement  results  from  two  contributions: 
one  from  the  "beam  deformation"  of  the  sample  and  one  from  the  local  deforma- 
tion of  the  mounting  pins.  The  relative  bias  due  to  sample  beam  deformation 
is  (Equation  (9a)): 


-(Ae  ) . 

I m z sample 

sample  e + (Ae  ) , 

o z sample 


Solving  for  (Ae  ) , and  substituting  e =0.1, 

3 z sample  3 o 


r 1 ,1 

r 1 ,] 

eRE^sample 

/ *RE^ sample 

(Ae  ) . 

z sample 


As  discussed  in  3.1.1,  this  strain  bias  corresponds  to  a force  Pq  (the 

clip  gage  reaction  force  at  10  percent  gage  strain).  For  the  same  force 

P , the  strain  bias  due  to  pin  deformation  is: 
o’ 


(Ae  ) . 
z pin 


••  M ...  ■ 
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The  net  strain  bias  is  therefore 


(Ae  ) 
z net 


(Ae  ) . + (Ae  ) . 

z sample  z pin 


(73) 


and  the  strain  experienced  by  the  sample  when  the  gage  strain  is 


e = 0. 1 is 
o 


e . = e + (Ae  ) .. 

sample  o z net 


Finally,  the  net  relative  bias  in  the  strain  gage  measurement  is 


^eRE^  net 


e - e . 
o sample 

£sample 


0.1 


'sample 


- 1 


(74) 


To  calculate  the  combined  measurement  bias,  the  calculator  program  used  to 

generate  the  results  in  3*1  was  modified  slightly  to  allow  calculation  of 

(e„r)  . for  fixed  values  of  P and  varying  modulus;  the  calculator  pro- 

Rt  sample  o 

gram  for  the  pin  deflection  analysis  was  used  to  calculate  the  term  (Ae  /P)  . 

z p i n 

for  the  same  modulus;  and  a new  calculator  program  was  written  to  carry 
out  the  calculations  of  Equations  ( 7 1 ) » (72),  (73),  and  (.lb). 


The  results  for  two  selected  values  of  PQ  are  shown  in  Figures  2b  and  25. 
Each  set  of  computations  was  carried  out  for  an  eccentricity  of  zero  and 
an  eccentricity  of  0.050  inches. 


The  value  of  PQ  used  in  Figure  2b  was  0.322  lb  (corresponding  to  a spring 
constant  of  k = 8.05  Ib/in  at  the  fixed  gage  length  of  0.4  inches).  This 
corresponds  to  the  calculated  stiffness  of  Thiokol's  early  experimental 
clip  gage,  which  was  residual  hardware  from  the  Simulated  Air  Launch  En- 
vironment (SALE)  program  (Reference  7).  This  gage  is  the  "baseline"  gage 
in  3.4.  It  appears  that  this  gage  would  produce  a measurement  bias  that 
would  be  unacceptable  for  most  engineering  purposes  over  the  entire  range 
of  propellant  modulus  considered.'  Since  the  gage  was  considered  too  stiff 


1.  Thiokol  reported  (Reference  6)  a negative  bias  of  approximately  35  percent 
for  the  baseline  gage  used  on  a propellant  sample  at  room  temperature  in  low 
rate  tests.  This  is  consistent  with  Figure  24. 
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U BIAS  in  measured  strain) 


SAMPLE  BENDING  ALONE 


SAMPLE  BENDING  PLUS  PIN  DEFLECTION 


^ SAMPLE  BENDING  PLUS  PIN  DEFLECTION 
(e=0.050  in) 


10  '“PROPELLANT  ’ ' 


Figure  24.  Expected  Strain  Measurement  Bias  for  Baseline 
Gage,  Used  on  Tensile  Specimen  of  Figure  2 
(L  =0.40  in) . 
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Figure  25 


Expected  Strain* Measurement  Bias  for  a 
Gage  With  a Spring  Constant  of  2.56  lb/ 
Used  on  Tensile  Specimen  of  Figure  2 (L 


4 
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for  the  intended  use,  Thiokoi  attempted  to  develop  a much  more  compliant 
gage.  The  final  prototype  gage  produced  had  a stiffness  of  2.56  lb/in 
and  a gage  length  of  0.40  inches.  Figure  25  shows  the  predicted  measurement 
bias  for  this  gage.  Comparison  of  Figures  2**  and  25  reveals  that  the  68 
percent  reduction  in  gage  stiffness  reduces  the  measurement  bias  by  a 
factor  of  roughly  50  percent  at  100  psi  and  somewhat  more  at  the  higher 
values  of  propellant  modulus.  The  remaining  measurement  bias,  however, 
is  still  much  larger  than  would  be  acceptable  for  most  engineering  purposes. 


Comparison  With  Finite  Element  Analysis  Results. 


As  a check  on  the  validity  of  the  models  used  in  the  preceding  analyses, 
independent  finite  element  analyses  of  the  mounting  pin/propellant  sample 
system  were  made  with  the  TEXGAP-30  structural  analysis  code.  The  finite 
element  model  used  the  sample  geometry  of  Figure  1,  with  Lq  » 0.40  inches 
and  a clip  gage  eccentricity  of  0.050  inches.  The  pin  diameter  and  material 
properties  were  the  same  as  used  in  this  report.  Three  values  of  propellant 
modulus  were  used:  200  psi,  1000  psi,  and  20,000  psi. 

The  deflection  of  the  pin  for  the  propellant  modulus  of  1000  psi,  as  pre- 
dicted by  the  finite  element  analyses  and  by  the  models  used  in  this  report, 
is  shown  in  Figure  26.  The  total  deflection  predicted  by  the  latter  is 
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Modifying  Equation  (7)  to  account  for  the  loading  eccentricity  e, 


I.  These  analyses  were  performed  by  Capt.  J.  Hildreth  in  the  process  of 
checking  out  the  new  code  on  AFRPL's  COC  6400  computer  system. 
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Figure  26.  Predicted  Total  Deflection  of  Mounting  Pins 
Including  Both  Sample  Bending  and  Local  De- 
flection (Propellant  Modulus  * 1000  psi). 
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The  "beam"  strain  at  the  sample  surface  nearest  the  clip  gage  (point  A in 
Figure  26)  Is 
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Similarly,  the  "beam"  strain  at  point  B is 
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and  the  sample  beam  deformation  contribution  to  the  "strain"  at  the  gage 
location  (point  A1  in  Figure  26)  is 
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These  relationships,  evaluated  at  P * I lb,  £ = 1,000  psi,  and  e ” 0.050  inches, 
result  in  the  deflection  labeled  "Sample  Flexure"  in  Figure  26.  To  obtain  the 
total  deflection  curves  labeled  "Rigid  Beam",  "Infinite  Beam",  and  "Interpolated 
Result"  in  Figure  26,  the  local  component  of  deflection  for  each  case  (shown 
in  Figure  22)  was  added  to  the  "Sample  Flexure"  deflection  at  each  point. 

The  results  show  excellent  agreement  with  the  TEXGAP-30  analysis  results, 
as  can  be  seen  in  Figure  26.  At  the  clip  gage  location,  the  two  results  agree 
within  2.5  percent. 

A further  comparison  with  the  TEXGAP-3D  results  is  given  in  Figure  27.  The 
deflections  from  the  finite  element  analyses  were  converted  to  a relative  bias 
in  strain  measurement  by  the  equation 
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Us i ng  Pq  ■ 0.32219  lb  and  Lq  « 0.4  lb,  the  three  points  plotted  In  Figure 
27  were  obtained,  showing  satisfactory  agreement  with  the  results  obtained 
from  the  models  used  in  this  report  over  the  wide  modulus  range  considered. 
It  appears  likely  that  correction  of  the  slight  discrepancy  in  calculating 
the  propellant  "foundation  stiffness"  (see  the  footnote  to  the  discussion 
of  Figure  20)  would  improve  the  agreement  between  the  models  in  this  report 
and  the  TEXGAP-3D  analysis. 
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3.  Clip  Gage  Structural  Behavior. 

The  analyses  in  this  subsection  were  conducted  to  estimate  the  stiffness 
and  relative  sensitivity  of  clip  strain  gages  and  then  to  use  these  results, 
in  combination  with  those  of  3.1,  to  explore  trade-offs  between  measurement 
bias  and  gage  output.  As  discussed  at  the  beginning  of  section  3,  the 
analyses  of  the  present  subsection  were  completed  before  the  importance  of 
the  local  deflections  near  the  mounting  pins  was  realized  and  before  the 
analyses  of  3-2  and  3.3  were  performed.  Except  where  otherwise  noted, 
therefore,  the  strain  measurement  bias  as  used  in  this  subsection  includes 
the  sample  beam  deformation  component  only.  Based  on  the  results  of  3-3, 
the  measurement  bias  actually  expected  would  be  roughly  twice  the  bias 
based  on  sample  beam  deformation  alone. 

3.^.1  Structural  Model. 

The  deflection  analysis  of  the  clip  gage  structure  was  based  on  simple 
elastic  beam  theory.  In  contrast  to  the  propellant  sample  analyzed  in 
the  preceding  subsections,  the  basic  structure  of  the  clip  gage  (see 
Figure  28)  is  composed  of  a material  which  can  be  regarded  as  being 
linearly  elastic  with  practically  no  error.  The  gage  structure  is  fabri- 
cated from  Beryllium  copper  foil  with  Young's  modulus  of  18.5  X 106  psi 
(Reference  8).  The  strain  sensors  themselves  were  not  modeled  in  the  present 
analysis;  the  actual  gages  can  therefore  be  expected  to  be  somewhat  stiffer 
than  predicted.  The  clip  gage  geometry  considered  is  shown  in  Figure  6. 
Several  different  mounting  techniques  have  been  used  to  attach  clip  strain 
gages  to  free  propellant  surfaces,  including  direct  bonding  to  the  pro- 
pellant surface  with  an  adhesive,  suspension  upon  pins  attached  in  various 
ways  to  the  propellant  surface  (as  in  Figure  2 and  Figure  7),  and  combinations 
of  pins  and  adhesive  bonds.  For  the  purpose  of  this  analysis,  only  the 
mounting  techniques  similar  to  those  of  Figures  2 and  7 were  considered. 

The  principal  reason  for  this  restriction  is  that  other  common  mounting 
techniques,  which  constrain  the  "feet"  of  the  strain  gage  to  zero  rotation, 
greatly  increase  the  stiffness  of  the  gage. 
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Figure  28.  Clip  Gage  Geometry  Used  in  Deflection  Analysis 


To  analyze  the  gage  structure,  it  was  divided  into  the  four  free  bodies 
shown  in  Figure  29.  The  deflection  equations  for  each  free  body  were 
written  out  and  the  resulting  total  deflection  and  rotation  at  points  A, 

B,  C,  and  D were  found  by  superposition.  The  actual  calculations  were 
performed  using  a calculator  program  written  with  subroutines  for  straight 
and  curved  beams  and  the  appropriate  logic  for  calculating  beam  dimensions 
and  applied  loads  from  the  overall  gage  geometry.  The  basic  calculation 
determines  gage  deflection  for  the  unit  load  applied  at  point  0 (Figure 
29).  A subsequent  calculation  determines  the  force  P^  required  to  obtain 
a strain  of  10  percent  over  the  gage  length  L . 

The  straight  beam  sections  of  Figure  29  are  each  treated  as  cantilevered 
beams  with  no  correction  for  "wide  beam"  Poisson  effect  (a  check  was  made 
to  verify  that  this  effect  was  negligible).  The  curved  beams  are  treated 
as  curved  cantilevered  beams  with  the  neutral  surface  halfway  between  the 
inner  and  outer  radii  of  the  beam;  the  Poisson  effect  was  again  ignored. 
(These  modeling  assumptions  probably  cause  noticeable  error  in  the  curved 
beam  deflection  itself;  however,  the  deflections  contributed  by  the  curved 
beams  are  small  relative  to  the  total  deflection,  and  the  relative  error  in 
the  total  deflection  is  assumed  to  be  negligible.)  The  "fixed"  support  for 
each  cantilever  beam  is  located  at  the  end  nearest  the  origin  of  coordinates 
in  Figure  29.  The  deflection  calculations  are  carried  out  first  for  the 
horizontal  beam  anchored  at  the  origin  of  coordinates  (Figure  29),  producing 
the  calculated  vertical,  horizontal,  and  rotational  deflections  at  point  A. 
These  deflections  are  then  considered  as  a motion  of  the  "fixed"  support  of 
the  curved  beam  AB;  the  rigid-body  motion  of  the  curved  beam  corresponding 
to  the  motion  of  the  support  plus  the  deflections  of  the  curved  beam  are 
combined  to  determine  the  vertical,  horizontal,  and  rotational  motion  of 
point  B.  This  procedure  is  continued  until  the  motion  of  point  D is  deter- 
mined. 


; 
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3.4.2  Analysis  Results. 

The  results  of  the  clip  gage  analyses  using  the  model  described  above  are 
discussed  In  the  following  subsections. 

3.4.2. I Effects  of  Gage  Geometry  on  Stiffness  and  Sensitivity. 

Figure  30  shows  the  variation  of  the  parameter  Pq  (force  at  10  percent  gage 
strain)  with  clip  gage  height.  Figure  31  shows  the  effect  of  changes  in 
clip  gage  height  on  a measure  of  the  gage  output  called  the  relative  sensi- 
tivity, defined  as  the  ratio  of  strain  in  the  clip  gage  center  section  to 
strain  in  the  sample  gage  length  Lq.  Figures  30  and  31  show  that  both  the 
gage  stiffness  and  relative  output  increase  rapidly  as  the  height  decreases. 
These  figures  also  show  that  changes  in  radius  of  the  corner  bends  in  the 
clip  gage  structure  have  only  a minor  impact  on  stiffness  and  output.  It 
is  also  apparent  from  the  two  figures  that  changes  in  height  and  bend  radius 
tend  to  preserve  a direct  relationship  between  output  and  gage  stiffness: 
a decrease  in  stiffness  is  obtained  at  the  expense  of  a reduction  in  gage 
output  for  a given  strain.  Finally,  comparison  of  Figures  30  and  31  shows 
that  a reduction  in  the  thickness  of  the  clip  gage  from  0.006  inches  to  0.004 
inches  is  predicted  to  change  the  stiffness  much  more  than  it  changes  the 
output  (this  is  a highly  desirable  result).  These  qualitative  observations 
provided  the  incentive  to  investigate  the  trade-offs  between  gage  sensitivity 
and  measurement  bias  (produced  by  the  clip  gage  stiffness)  which  can  be 
realized  through  changes  in  clip  gage  configuration. 


Before  examining  output/accuracy  tradeoffs,  however,  the  specific  results 
shown  in  Figures  30  and  31  for  the  "baseline  configuration"  will  be  dis- 
cussed. The  baseline  configuration  corresponds  to  Thiokol's  Prototype 
Gage  No.  1 and  Prototype  Gage  No.  2 (Reference  9) . These  two  prototype 
gages  were  actually  successive  modifications  of  a clip  strain  gage  provided 
by  Konigsberg  Instruments,  Inc.  for  use  in  the  Simulated  Air  Launch  Environ- 
ment (SALE)  program  (AFRPL  Contract  F0461 1-73_C-0015) - This  particular 
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clip  gage  was  not  used  in  the  SALE  program  because  of  damage  which  rendered 
the  semiconductor  strain  gage  bridge  inoperative.  Although  the  gage  was 
nonfunctional,  Thiokol  considered  it  to  be  a useful  point  of  departure  for 
their  own  prototype  gage  development  (Reference  8).  Since  Thiokol  preferred 
to  use  foil  strain  gages  rather  than  semiconductor  strain  gages,  the  original 
semiconductor  gages  were  removed,  and  the  gage  structure  was  instrumented 
with  a four-active-arm  bridge  foil  strain  gage  arrangement  (Reference  8), 
resulting  in  Prototype  Gage  No.  1.  When  problems  with  the  strain  gage 
bonds  were  encountered,  the  foil  strain  gages  were  stripped  off  and  the 
gage  structure  was  re- instrumented  (again)  with  foil  strain  gages,  result- 
ing in  Prototype  Gage  No.  2 (Reference  9). 

Thiokol  measured  the  stiffness  of  this  gage  structure  in  the  various  stages 
of  modification  described  above;  the  results,  in  terms  of  a computed  gage 
stiffness  or  spring  constant,  are  shown  in  Table  1.  Table  I also  shows 
stiffnesses  calculated  from  Figure  30  and  from  Thiokol 's  own  independent 
beam-theory  analysis.  The  rather  large  differences  in  experimental  stiff- 
nesses between  Prototype  No.  2 and  Prototype  No.  1 (the  latter  having  the 
same  stiffness  as  the  Konrgsberg  gage)  has  not  been  explained.  It  is  evident, 
however,  that  the  model  used  to  produce  Figure  30  agrees  closely  with  Thiokol' 
analysis,  does  reasonably  well  at  predicting  the  behavior  of  the  bare  gage 
structure  (within  about  10  percent),  and  probably  substantially  under-pred icts 
the  stiffness  of  a "real"  clip  gage.  The  latter  observation  is  to  be  expected 
since  the  stiffness  added  by  the  foil  (or  semiconductor)  strain  gages,  wires, 
and  epoxy  adhesive  is  not  included  in  the  structural  model.  The  model  does 
appear  to  be  sufficiently  accurate  for  the  gage  design  trade-off  analyses 
of  the  following  subsection. 

3 . A . 2 . 2 Sensitivity  Versus  Measurement  Bias. 

The  incentive  for  seeking  improvement  in  the  measurement  bias  is  readily 
apparent  from  Figure  24  (3.3.1),  which  predicts  a very  large  negative  bias, 
even  if  the  effect  of  local  deflections  of  the  mounting  pins  is  ignored. 

Since  Thiokol  was  using  metal  foil  strain  gages  (which  have  a small  gage 


TABLE  1.  STIFFNESS  OF  EARLY 
THIOKOL  PROTOTYPE  GAGES  (MODIFIED  KONIGSBERG  GAGE) 


NOMENCLATURE 

SOURCE  OF  DATA 

CALCULATED  STIFFNESS 

Original  Gage 

Thiokol -Measured  force  at 
0.030  inches  extension 
[Reference  8] 

10.3 

1 b/  i n 

Gage  Structure  Only 

Thiokol -Measured  force  at 
0.030  inches  extension 
[Reference  8] 

7.3 

lb/  i n 

Gage  Structure  Only 

Thiokol-Analytical  predic- 
tion from  beam  theory 
[Reference  8]* 

8.0 

1 b/  in 

Gage  Structure  Only 

Figure  30* 

8.05 

1 b/  i n 

Prototype  Gage  No.  1 

Thiokol -Measured  force  at 
0.030  inches  extension 
[Reference  6] 

10.3 

1 b/  i n 

Prototype  Gage  No.  2 

Th iokol -Measured  force  at 
0.030  inches  extension 
[Reference  9] 

8.0 

lb/ in 

*Based  on  gage  length  of  0.400  inches 


factor  compared  to  semiconductor  strain  gages),  it  was  apparent  that  the 
sacrifice  in  output  to  achieve  decreased  gage  stiffness  (hence  decreased 
measurement  bias)  must  be  minimized.  The  trade-off  studies  discussed  below 
were  therefore  undertaken  to  acquire  the  insight  needed  to  assess  and  guide 
Thiokoi's  attempts  to  develop  clip  strain  gages  and  mounting  techniques 
which  would  allow  measurements  of  acceptable  accuracy  for  the  purposes  of 
the  Improved  Properties  program. 

Viewed  in  an  overall  sense,  the  problem  is  simply  posed.  The  strain  at  the 
strain  gage  location  must  be  relatively  increased  while  the  overall  stiffness 
of  the  clip  gage  is  relatively  decreased.  Optimally,  this  would  be  accom- 
plished by  making  all  parts  of  the  transducer  except  the  measurement  section 
as  rigid  as  possible  and  by  maximizing  the  strain  in  the  measurement  section 
for  a given  bending  stiffness  of  the  measurement  section.  Given  a design 
concept  that  accotnplsihes  these  objectives,  the  bending  stiffness  of  the 
gage  could  then  be  reduced  to  the  desired  level  (or  to  the  level  which  pro- 
duces the  minimum  acceptable  electrical  output).  A few  additional  practical 
considerations  would,  of  course,  be  present  (e.g.,  elastic  instability, 
gage  durability  in  routine  handling  and  use,  and  vibration  character istics 
under  dynamic  loads). 

Rather  than  attempting  to  provide  this  "ultimate  solution",  however,  it 
was  considered  more  practical  to  observe  the  limitations  imposed  by  Thiokol 
on  their  experimental  development.  These  were,  In  essence: 

(a)  Retain  the  same  (gage  structure)  material; 

(b)  Retain  the  foil  strain  gages  (as  opposed  to  semiconductor 
gages) ; 

(c)  Keep  the  gage  "simple  to  manufacture." 

Three  potential  techniques  for  reducing  the  gage  measurement  bias  are 
explored  in  Figure  32:  reducing  the  gage  thickness,  increasing  the  gage 
height,  and  decreasing  the  gage  length  (Lq  in  Figure  l).  The  two  curves  in 
Figure  32  labeled  t * O.OQl*  and  t = 0.006  were  obtained  by  cross  plotting 
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the  relative  sensitivity  for  a given  value  of  H (from  Figure  31)  against 
the  relative  bias  (from  Figure  3)  corresponding  to  the  value  of  PQ  produced 
by  that  value  of  H (from  Figure  30).  The  dashed  line  representing  the 
effect  of  a reduction  in  gage  length  was  obtained  in  a similar  manner,  as 
shown  in  Figure  33.  Note  that  the  force  for  a -1  percent  bias  is  a function 
of  the  gage  length  (the  data  in  Figure  33  are  obtained  from  Figure  4).  The 
force  required  for  a -I  percent  bias  was  found  to  intersect  the  curve  of 
Pq  versus  gage  length  at  a gage  length  of  approximately  0.2  inches;  this 
value  of  the  gage  length  produces  a relative  output  of  0.059.  This  value 
is  plotted  in  Figure  32  as  the  end  point  of  the  dashed  line.  (The  gage 
height  of  0.25 0 inches  was  selected  arbitrarily). 

It  is  evident  from  Figure  32  that  reducing  the  gage  thickness  results  in 
less  penalty  (in  terms  of  reduced  output)  for  the  resulting  decrease  in 
measurement  bias  than  does  a gage  height  increase.  In  turn,  increasing 
the  gage  height  produces  less  penalty  for  the  measurement  bias  reduction 
than  does  reducing  the  gage  length.  Figure  32  also  shows  that  none  of  the 
three  techniques  considered  allows  a reduction  in  the  measurement  bias 
without  a sacrifice  of  output  level. 


It  would  appear,  as  discussed  above,  that  a significant  improvement  in 
the  sensitivity-versus-bias  character i s t i cs  of  the  clip  gage  might  be 
obtained  if  the  gage  legs  could  be  made  rigid.  For  a gage  with  rigid 
legs,  the  deformation  would  all  occur  within  the  horizontal  section,  as 
illustrated  in  Figure  3Mb).  From  Figure  3Ma)  , the  constant  bending 
moment  in  this  section  is  seen  to  be  PQh.  The  relative  sensitivity  is 


Mc/EI 

£ 

O 


P hc/EI 
o 


From  Figure  3Mb),  the  strain  (for  small  values  of  0)  in  the  gage  length 

L i s 
o 
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Equations  (75)  and  (76)  were  used'  to  calculate  the  relative  sensitivity- 
versus-bias  curve  for  a rigid-leg  gage  having  the  dimensions  and  modulus 
of  the  baseline  gage.  The  result  is  compared  in  Figure  35  with  the  sensi- 
tivity-bias characteristics  of  the  baseline  gage.  As  expected,  the  rigid- 
leg  gage  provides  a higher  relative  sensitivity  for  a given  measurement 
bias  but  is  stiffer  than  a flexible-leg  gage  of  the  same  dimensions. 

Thiokol  produced  a prototype  gage  which  approached  the  ideal  rigid-leg 
configuration.  This  gage,  which  was  the  final  prototype  discussed  in 
3.3-1,  had  a measured  spring  constant  of  2.56  Ib/in.  This  prototype  was 
designated  the  "stiff-leg  gage";  its  configuration  was  basically  as 
shown  in  Figure  36.  Calculations  show  that  the  behavior  of  the  stiff-leg 
gage  was  intermediate  between  a rigid-leg  gage  and  a flexible-leg  gage 
with  the  same  dimensions,  and  that  the  stiff-leg  gage  provided  a substantial 
reduction  in  measurement  oias  with  only  a small  reduction  in  sensitivity 
compared  to  the  baseline  gage.  However,  the  measurement  bias  of  nearly 
3 percent  (neglecting  bias  due  to  local  deformation,  which  will  raise  the 
total  bias  to  more  than  5 percent)  is  still  far  too  high.  (See  Figure  37.) 

Equations  (75)  and  (76)  were  used  to  calculate  trajectories  of  S versus 

A 

Pq  with  specific  parameters  allowed  to  vary  while  all  other  parameters 
were  held  fixed.  Examining  the  results,  shown  in  Figure  38,  demonstrates 
the  sensitivity  of  a rigid-leg  clip  gage's  output  (for  a given  value  of 
P ) to  the  gage  design  parameters.  It  is  apparent  in  Figure  38  that  vary- 
ing the  gage  width  (w)  or  the  elastic  modulus  (E)  is  the  most  efficient 
way  to  decrease  the  value  of  PD  (and  hence  the  gage  bias  due  to  gage- 
sample  interaction).  Note,  however,  that  interactive  effects  between 
the  design  parameters  could  provide  more  efficient  combined  parameter 

A A 

changes  to  reduce  P . (For  instance,  a decrease  in  w together  with  a 

O /s  A 

decrease  in  L would  reduce  P while  increasing  S.)  These  rigid-leg 

O /\  /s 

gage  results  provide  upper  bounds  on  S-versus-PQ  trends  for  flexible- 
leg  gages,  as  illustrated  by  the  solid  (triangle  and  square)  symbols  in 


1.  The  bias  was  determined  from  P using  the  techniques  described  in  3.3, 
but  the  local  mounting  pin  deformation  was  ignored. 
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Figure  38,  which  were  determined  from  Figures  30  and  31.  Note  that  the 
corresponding  rigid-leg  and  flexible-leg  points  lie  on  lines  of  constant 

A A 

S/F  ; this  is  always  true  for  gages  with  the  same  center  section  geometry 
o 

and  materials. 

Based  on  these  results  from  the  idealized  rigid-leg  gage  model,  a number 
of  potentially  achievable  changes  to  the  gage  structure  were  analyzed 
using  the  realistic  model  to  provide  an  idea  of  the  potential  for  improve- 
ment. The  results  are  presented  as  curves  A through  F in  Figure  39-  (In 
considering  the  impact  of  these  results,  recall  that  the  strain  measurement 
bias  shown  neglects  local  pin  deformation,  so  that  the  actual  total  bias 
would  be  about  twice  that  shown.  Also,  the  bias  at  100  ps i modulus  would 
be  nearly  five  times  the  bias  at  500  ps i modulus,  while  propellant  modulus 
values  greater  than  500  ps i would  reduce  the  bias  in  a roughly  proportional 
relationship).  Curve  A shows  the  effects  obtainable  by  increasing  the 
height  of  the  baseline  gage,  holding  all  other  parameters  constant.  Curve 
B shows  the  effects  obtainable  by  changing  the  gage  material  to  a titanium 
alloy  and  duplicating  the  geometries  considered  in  curve  A;  curves  C and 
D show  the  effects  obtainable  with  an  aluminum  alloy  and  a magnesium  alloy, 
considering  the  same  geometries  as  in  curve  A.  Note  (from  Table  2)  that 
the  range  of  moduli  available  in  this  group  of  materials  is  from  18.5  X 10^ 
psi  down  to  6.5  X 10^  psi.  Further  improvements  are  shown  when  the  foil 
thickness  is  reduced  to  . 00 4 inch  and  the  width  to  0.100  (curves  E for  the 
original  beryllium  bronze  material  and  F for  the  magnesium  alloy).  Finally, 
the  effects  of  changing  to  the  magnesium  alloy,  adopting  Thiokol's  "stiff-leg" 

configuration,  reducing  the  thickness  to  .004  inch  and  reducing  the  width 

2 

and  length  of  the  measuring  section  to  0.100  inch  are  shown  by  curve  G. 

These  changes  are  probably  achievable,  although  the  feasibility  of  using 
magnesium  alloy  for  the  gage  may  be  questionable  (titanium  and  aluminum 
alloys  have  been  used  in  propellant  stress  gages  and  strain  gages  respectively). 

The  conclusion  to  be  derived  from  Figure  39  is  that  it  is  potentially  possible 
to  produce  a clip  strain  gage,  using  foil  strain  gage  sensing  elements,  with 


2.  Calculations  for  curve  G used  the  rigid  leg  gage  model  with  a stiffness 
reduction  factor  (obtained  from  results  shown  in  Figure  37)  of  0.649. 


EFFECTS  OF  VARYING  HEIGHT  WITH 
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A TABLE  2 FOR  GAGE  DESCRIPTIONS) 
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TABLE  2.  GAGE  PARAMETERS  USED  IN 
CALCULATIONS  OF  FIGURE  39 


SYMBOL 

MATERIAL 

E,  psl* ** *** 

w.  i n 

t,  In  L, 

_in  Hir 

i n 

A 

BERYLLIUM 

BRONZE 

18.5XI06 

. 125 

.006 

.256 

.175 

B 

T 1 TAN  1 UM 

ALLOY  RC-130A 

16.0XI06 

.125 

.006 

.256 

.300 

C 

ALUMINUM 

ALLOY  24S-T6 

I0.5XI06 

• 125 

.006 

.256 

.288 

D 

MAGNESIUM 

ALLOY  AZ3IX 

6.5X106 

.125 

.006 

.256 

.212 

E 

BERYLLIUM 

BRONZE 

18. 5X1 06 

.100 

.004 

.256 

.181 

F 

MAGNESIUM 

ALLOY  AZ 3 1 X 

6.5X106 

.100 

.004 

.256 

. 102 

MAGNESIUM 

ALLOY  AZ3IX 

6. 5X1 06 

.100 

.004 

.100 

.161 

* Data  from  Reference  10. 

**  Calculated  result. 

***  "Stiff-leg"  gage  treated  as  rigid-leg  gage  with  stiffness  scaled  per 
experimental  stiffness  of  prototype  "stiff-leg"  gage. 
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an  acceptable  output  level  (comparable  to  that  of  the  Thiokol  stiff-leg 
prototype)  and  a bias  of  less  than  1 percent  at  a 500  psi  propellant  modulus, 
including  local  deformation  effects.  The  bias  at  100  psi  propellant  modulus 
would  be  at  a reasonable  level  of  about  5 percent.  The  primary  area  of 
doubt  would  be  the  usability  of  magnesium  alloy  or  another  low-modulus 
material  for  the  clip  gage  structure. 

Another  consideration  for  clip  gage  improvement  is  the  use  of  semiconductor 
strain  gages  instead  of  foil  strain  gages.  Semiconductor  gages  would  allow 
the  use  of  much  lower  values  of  relative  sensitivity  since  they  have  up  to 
75  times  higher  gage  factor  than  the  foil  gages.  However,  examination  of 
Figure  30  suggests  that  attempting  to  utilize  this  approach  to  reduce  the 
bias  by  increasing  the  height  of  a given  gage  structure  would  result  in 
rather  long  legs  before  much  advantage  is  gained.  This  is  borne  out  by 
calculations:  extending  curve  E (Figure  39)  to  a relative  sensitivity 
of  .0037  reduces  the  relative  bias  (including  local  deformation  effects) 
to  about  0.2  percent  at  500  psi  propellant  modulus  or  roughly  1.0  percent 
at  100  psi  propellant  modulus.  However,  the  required  leg  length  is  0.415 
inch  and  the  structure  would  be  somewhat  flimsy.  The  same  bias  level  with 
the  gage  design  approach  of  curve  G requires  a leg  length  of  0.510  inch. 

These  clip  gages  would  probably  be  within  the  limits  of  practicality. 


3.5 


Further  Analysis. 


The  analyses  in  this  subsection  were  conducted  to  investigate  the  effect  of 
clip  gage  prestrain  on  strain  measurement  bias;  to  incorporate  time-depend- 
ent propellant  modulus  Into  the  bias  calculations;  and,  finally,  to  compare 
the  time-dependent  results  with  Thiokol's  experimental  measurements  of  gage 
response  to  evaluate  the  theoret ical calculat Ions. 

3.5.1  Effect  of  Gage  Prestrain. 

In  all  earlier  analysis  in  this  repo-t,  the  clip  gages  were  assumed  to  be 
attached  to  the  sample  or  to  the  mounting  pins  in  such  a way  that  the  gage 
would  be  in  an  undeformed  state  when  the  propellant  sample  strain  over  the 
gage  length  was  zero.  That  is,  it  was  assumed  that  the  gage  was  mounted 
without  prestrain.  In  actual  use,  however,  the  clip  gages  are  almost 
always  mounted  in  a prestrained  condition  so  that  reversal  of  the  direction 
of  the  force  on  the  strain  gage  does  not  occur.  This  practice  avoids  certain 
mechanical  problems  in  the  use  of  the  gages. 


Prestrain  is  obtained  by  a mismatch  of  dimensions  between  the  gage  mounting 
points  (e.g.,  holes  in  the  gage  feet)  and  the  sample  mounting  points  (e.g., 
the  mounting  pins)  which  define  the  gage  length  for  the  strain  measurements. 
For  the  purposes  of  this  analysis,  assume  that  the  mounting  pins  are 
implanted  in  the  sample  a distance  Lq  (the  nominal  gage  length)  apart,  but 
that  the  mounting  holes  in  the  gage  feet  are  purposely  spaced  at  a distance 
Lq  - d apart.  The  gage  will  therefore  be  in  a prestrained  condition  when 
assembled  with  the  propellant  sample  and  mounting  pins  as  in  Figure  7.  If 
the  sample  and  pins  were  infinitely  stiff,  the  prestrain  would  be 


(e 


pre' nominal 


where  the  subscript  "pre"  indicates  "prestrain", 
are  not  infinitely  stiff,  both  the  sample  and  the 


Since  the  sample  and  pins 
gage  will  deform  by  an 
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amount  equal  to  (Ae  ) (L  ) because  of  the  force  required  to  prestrain  the 
pre  o 

gage.  This  prestrain  force  Is 


pre 


kg  [d  ‘ (Aepre)  Lo] 


where  k is  the  gage  stiffness  or  spring  constant.  However,  the  preload 
9 

force  is  also  given  by 


pre 


k (Ae  ) L . 
s pre  o 


Solving  for  Ppr£  and  epre  , 


pre 


k k d 

_J_2 

k + k 
s g 


Ae 


k d 
-1L 


pre 


(k  + k ) L 
s g o 


Recognizing  that  propellant  is  viscoelastic,  k$  will  be  different  for  a 
long-term  loading  situation  (prestrain)  than  for  a short  term  loading  of 
the  sample.  Define  the  stiffness  for  prestrain  as  (ks)pre-  Then 


pre 


(k  ) k d 

TrFrhr 


s pre 


and 


(Ae) 


k d 

_S 


pre 


[ (k  ) + k ] L 

s pre  g o 


The  experimental  practice  is  to  assign  a strain  measurement  of  zero  to  the 
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gage  output  at  this  condition.  Therefore  when  a strain  change  Ae  from 

gage 

this  condition  occurs,  the  measured  strain  value  is 


e = Ae  ■ e - e 

meas  gage  gage  gage 


pre 


gage 


[Lo  ' 


where  e is  the  actual  strain  (AL  /L  ) of  the  gage.  The  clip  gage 
gage  o o 

reaction  force  at  this  condition  is 


P - P + (AP)  - P + k (Ae  ) L 

pre  e pre  g meas  o 

meas 


The  total  "bias"  strain  produced  by  P is 


or 


Ae  * (Ae  ) + (Ae  ) • e . - e 

z z z sample  gage 

pre  meas  K 


Ae 


e . - Ae  + 

sample  gage 


[Lo  " (A£z)pre] 


so 


Ae 


meas 


Ae  . - Ae  = Ae  . - A(e  ) 

sample  gage  sample  meas 


since  the  term  in  brackets  cancels  out.  Substituting  for  A(e  ), 

meas 


Ae 


z k L 

meas  s o 


k (Ae  ) 
g meas 


or 


Ae 


meas 


k 

_a. 

K 

s 


e . - Ae 

sample  z 


meas 
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or 


Ae 

z 

meas 


* - £ 

k + k sample 

s g 


Now  the  gage  stiffness  in  the  terminology  of  3.1  is 


k 

9 


P 

o 


c L 


o o 


A P 

o 

Ae  L 


o o 


and  the  sample  stiffness  is  (for  the  geometry  of  Figure  1) 


I 4 . Vo  L 


so  the  relative  measurement  bias  becomes 


Rc 


meas 


-[A  Po  /(Aeo  Lo)] 
E A P 

1 4 . 4ol  ° 


-14.4  A P 

o 

Ae  E + A P 
o o 


Ae  L 
o o 


We  are  therefore  led  to  the  conclusion  that  prestrain  does  not  affect  the 
measurement  bias  in  the  usual  method  of  clip  gage  usage J 


3.5.2  Incorporation  of  Viscoelastic  Propellant  Behavior. 

All  previous  analyses  in  this  report  have  treated  the  propellant  sample  as 
a linear  elastic  body.  Time-dependent  effects  will  now  be  incorporated  into 
the  analytical  model  by  including  consideration  of  viscoelastic  propellant 
response. 

Thiokol's  characterization  of  the  TP-H8260  propellant  used  in  clip  gage 
experiments  under  the  Improved  Properties  program  provided  the  relaxation 
modulus  data  shown  in  Figure  40.  The  calculator  programs  discussed  in 


1.  However,  a viscoelastic  analysis  shows  short  term  effects  (see  Figure  49). 
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3.3.1  were  used  to  calculate  the  strain  measurement  bias  for  the  Thiokol 
stiff-leg  gage  (described  in  Figure  36)  at  selected  time  intervals,  using 
the  stress  relaxation  modulus  given  at  each  time  by  Figure  40.  These 
calculations  included  both  the  bias  due  to  sample  beam  deformation  and 
the  bias  due  to  local  deformation  of  the  mounting  pins.  Calculations  were 
carried  out  for  two  loading  eccentricities:  e ■ 0 and  e = 0.050  in. 

The  results,  shown  in  Figure  1*1,  give  the  time-dependent  response  of  the 
clip  gage  (due  to  gage-sample  interaction)  to  a step  function  strain  applied 
to  the  sample  at  t/a^  ■ 0,  where  t is  time  and  a^.  is  the  t ime-temperature 
(WLF)  shift  factor.  The  response  is  shown  in  terms  of  the  relative  strain 
measurement  bias. 


As  shown  in  Figure  41,  the  time-dependent  relative  measurement  bias  is 

closely  approximated  by  a linear  function  of  Log  t/aT  over  the  range 

”2  2 

of  interest  for  present  purposes  (10  t/a^  <_  10  ).  Having  this 
approximate  relationship,  i.e.  (taking  a^.  = 1) 


-eRr  = A + B Log  t. 


the  response  at  time  t to  a step  strain  input  Ae . applied  at  time  (t  ).  is 

• o i 

(A  Response).  ■*  (I  + e^)  Ae(  » (1-A-B  Log  t)  AEj. 

Therefore,  for  a strain  history  made  up  of  a series  of  superposed  step 
functions  in  time, 


n n T 1 

5.  A Response  - Z 1-A-B  Log  (t-t  ).  Ac,  6, 

i«l  i-1  L o .J  i , 


(78) 


where 


A . - I for  (t-t  ) , > 1 

i o I — 


0 for  (t-t  ).  < I 
O I 
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A calculator  program  was  written  to  carry  out  the  calculation  of  Equation 
(78)  for  a series  of  five-step  functions  of  arbitrary  levels  occurring  at 
equal  t i me  intervals  (response  Is  output  at  time  increments  of  (t-tQ)/IO; 

(t-t  ) is  assumed  to  he  equal  to  five  units  of  time). 

Using  the  calculator  program,  tho  response  to  the  applied  strain  history 
is  shown  in  Figure  A2.  This  strain  history  was  selected  to  agree  with  that 
used  by  Thiokol  in  a laboratory  gage  calibration  test  on  a gaged  uniaxial 
propellant  sample  (Reference  12).  The  strain  increments  correspond  to 
Thiokol's  reported  data;  the  time  interval  between  steps  (5  minutes)  and 
the  interval  between  strain  application  and  gage  data  acquisition  (5  minutes) 
were  selected  based  on  the  best  estimate  of  the  test  personnel,  since  this 
data  was  not  recorded. 

3 . 5 . 3 Compar i son  With  Test  Results. 

Ihe  lhiokol  laboratory  test  of  the  gaged  sample  (Reference  12)  was  carried 
out  with  considerable  attention  to  detail.  The  "true"  sample  strain  at 
each  strain  step  was  determined  by  travelling-microscope  observation  of 
the  mounting  pins  at  the  gage  attachment  points  with  the  qage  removed. 

The  sample  was  allowed  to  "rest"  at  zero  strain,  then  the  gage  was  mounted 
on  the  pins  (as  shown  in  Figure  7)  and  the  strain  history  was  repeated. 

During  the  second  strain  history,  both  the  electrical  strain  gage  output 
and  the  pin-to-pin  optical  measurements  were  recorded.  The  results,  plotted 
in  Figures  V}  and  A4  (solid  symbols),  were  determined  by  averaging  five 
individual  observations  of  each  data  point. 

Interestingly,  two  factors  which  are  observable  in  Figure  ^3  produced  a 
false  sense  of  security  in  earlier,  less  thorough  tests.  The  first  is  the 
increased  strain  in  the  gage  length  of  the  sample;  the  second  is  the 
discrepancy  between  the  gage  electrical  output  and  the  actual  strain 
imposed  on  the  gage.  Both  factors  produced  errors  which  compensated  in 
part  for  the  inherent  bias  due  to  the  gage/sample  interaction. 
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Figure  44.  Comparison  of  Thiokol  Calibration  Experiment  Data 
With  Analysis  Predictions  (Driving  Strain  - Gaqe 
Length  Strain  with  Nb  Gage)  (Ref.  12). 
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(The  error  introduced  by  the  increased  strain  in  the  gage  length  was  present 
because  this  data  was  not  observed  in  the  earlier  tests.)  To  Illustrate, 
ignoring  these  two  factors  would  lead  to  an  estimate  of  the  measurement  bias, 
at  5.6  percent  applied  strain,  of  about  -7  percent  of  the  "true"  value. 
Considering  the  two  factors  leads  to  the  conclusions  that  (1)  the  in-situ 
calibration  of  the  gage  does  not  agree  with  Its  original  calibration;  and  (2) 
in  the  absence  of  this  compensating  error,  the  measurement  bias  at  5.6  per- 
cent applied  strain  would  be  about  ”23  percent  of  the  "true"  value. 

The  analysis  results  of  Figure  42  (e  ■ 0.050  in)  and  similar  results  for 
e » 0 are  plotted  in  Figure  1*3-  The  results  for  e **  0 are  quite  close  to 
the  optically  observed  results,  while  the  results  for  e ■ 0.050  inch  fall 
below  the  observed  results.  It  was  noted,  however,  that  the  higher  strain 
near  the  center  of  the  sample  length  is  not  accounted  for  in  the  gage 
analysis.  The  analysis  was  repeated  using  the  observed  (no  gage)  strain 
in  the  sample  gage  length  as  the  "driving"  strain.  The  results,  shown 
in  Figure  44,  indicate  that  the  value  of  e - 0.050  inch  produces  better  agree- 
ment with  observed  data  than  does  e ■ 0.  In  fact,  e » 0.050  inch  is  the 
"best  guess"  of  the  actual  loading  eccentricity. 

In  all,  the  very  close  agreement  between  analytical  and  experimental  results 
shown  in  Figure  44  provide  considerable  confidence  in  the  analyses  and 
thence  in  the  conclusions  reached. 

3.5.4  Effects  of  Cyclic  Loads  And  Prestrain. 

The  calculator  program  used  for  the  analysis  discussed  in  the  preceding 
paragraphs  was  used  to  investigate  the  effects  of  cyclic  loads  in  a very 
limited  manner.  Three  different  strain  histories,  all  Involving  cycling 
about  a mean  value  in  a stepwise  fashion,  were  applied  to  the  model;  the 
r-aults  are  shown  in  Figures  45  to  48.  In  Figure  46,  note  that  the  hystere- 
sis is  reduced  by  selecting  a longer  interval  between  application  of  strain 
and  data  acquisition;  however,  this  does  not  always  reduce  the  measurement 
bias. 
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Figure  48.  Apparent  Strain 


Ratcheting 


Propellant  Response  (Strain  History  No.  3) 


An  even  more  interesting  effect  is  seen  in  figure  A8.  This  is  a strain 
"marching"  or  "ratcheting"  effect,  which  closely  resembles  results  obtained 
in  certain  instrumented  motor  tests  (References  13  and  14).  Whether  the 
viscoelastic  effects  of  gage/propellant  interaction  seen  in  Figure  A8  has 
any  bearing  on  the  instrumented  motor  results,  which  were  interpreted  as 
real  phenomena,  is  unkown.  However,  some  investigation  of  this  point  may 
be  worthwhi le. 

To  investigate  the  effects  of  prestrain  on  gage  measurement  bias,  the 
following  scenario  was  utilized: 


At  time  tj  - tQ,  a step  prestrain  (0.01  in/in  was  used  in 
the  analysis)  is  applied  to  the  sample. 


At  time  tj,  a step  strain  of  magnitude  Ae  is  applied  in 
the  same  direction  as  the  prestrain. 


At  time  t (in  the  analysis,  t = tj  + 5),  the  strain  sensed 
by  the  gage  is  recorded.  The  strain  sensed  by  the  gage  at 
time  tj  (i.e.,  just  before  application  of  the  strain  Ae)  is 
subtracted  from  the  recorded  value  to  obtain  the  measured 
strain  corresponding  to  Ae.  For  a linear  viscoelastic 
sample  material,  this  is  entirely  equivalent  to  the  clip 
gage  usage  described  in  3.5.1. 


In  order  to  carry  out  the  calculations  for  prestrain  times  greater  than 
20  minutes,  the  time  variable  t in  the  calculator  program  described  in 
3.5.3  was  replaced  by  0 * (t/n)  which  required  that  the  value  of  A (see 
Equation  (78))  be  replaced  by  A + B Log  n and  that  the  time  value  output  be 
multiplied  by  n. 
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The  results  of  the  prestrain  analysis  are  shown  In  Figure  49.  Evidently 
a relaxation  time  of  an  hour  or  so  is  ample  for  all  prestrain  effects  to 
disappear  in  gage  response  at  5 minutes  after  loading.  C .culations  for 
gage  response  at  0.1  minutes  after  loading  show  negligible  prestrain 
effects  even  after  only  5 minutes  of  relaxation. 
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CONCLUSIONS  AND  RECOMMENDATIONS. 


A. 


A. I Genera)  Conclusions . 


This  report  covers  extensive  analyses  of  the  effects  of  gage/propellant 
interaction  on  the  accuracy  of  propellant  strain  measurement  using  clip 
strain  gages.  The  analyses  directly  address  strain  measurement  on  uniaxial 
tensile  specimens  of  a particular  configuration  used  in  the  AFRPL/Th iokol 
Improved  Properties  program.  However,  as  discussed  below,  the  conclusions 
are  true  for  other  applications  as  well. 

The  following  general  conclusions  are  supported  by  the  analysis  results: 

a.  At  least  some  state-of-the-art  clip  strain  gages  are  too 
stiff  for  accurate  measurement  of  solid  propellant  strain. 
The  starting  point  of  Thiokol's  gage  development  efforts 
was  a clip  gage  structure  supplied  by  Konigsberg  Instruments 
and  originally  intended  for  use  on  the  AFRPL/Thiokol  SALE 
program.  This  structure  was  stiff  enough  to  introduce 
enormous  measurement  bias.  Even  Thiokol's  best  effort  at 
producing  a compliant  gage  induced  a strain  bias  of  more 
than  20  percent  of  the  measured  strain. 


b.  Drastic  changes  in  clip  gage  materials  and  design  could 
potentially  eliminate  or  reduce  measurement  bias  due  to 
gage/propellant  interaction. 

c.  The  strain  gage/propellant  interaction  problem  can  be 
analyzed  with  acceptable  accuracy.  The  model  used  in  the 
analyses  in  this  report  agrees  extremely  well  with  both 
3-D  finite  element  analyses  and  with  experimental  measure- 
ments. 
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As  a result  of  early  predictions  of  the  analyses  discussed  in  this  report 

and  Thiokol's  confirming  laboratory  tests  of  gaged  propellant  samples,  a 

decision  was  made  to  abandon  clip  strain  gages  as  a measurement  tool  in 

that  program.  Alternate  techniques  were  developed  for  measurements 

allowing  the  determination  of  propellant  Poisson's  ratio  (these  techniques 

are  discussed  in  Reference  1).  i 

4 . 3 Implications  For  Other  Applications. 

A careful  examination  of  the  cons iderat ions  involved  (see  3-2.4. 4)  in 
combining  the  bias  due  to  local  mounting  pin  deformation  with  that  due  to 
sample  "beam  deformation"  leads  to  the  (perhaps  surprising)  conclusion 
that  the  gage-propellant  interaction  effects  are  very  nearly  the  same  for 
a clip  gage  mounted  on  a semi  - inf  in i te  slab  as  for  a clip  gage  mounted  on  a 
uniaxial  sample.  This  conclusion  applies  equally  well  to  a clip  gage 
mounted  on  the  bore  surface  of  an  instrumented  motor.  It  is  clear,  there- 
fore, that  action  is  needed  to  provide  acceptable  means  of  experimental 
propellant  strain  measurement. 

4.4  Recommendations . 

Experimental  measurements  of  stress  and  strain  are  necessary  for  understanding 
of  propellant  mechanical  behavior,  both  in  simple  laboratory  samples  and  in 
actual  rocket  motors.  The  tools  and  techniques  for  making  these  measurements 
must  be  developed,  or  the  lack  of  such  aids  will  eventually  become  an  ob- 
stacle to  further  improvements  in  theoretical  modeling  of  propellant  be- 
havior as  well  as  experimental  investigation.  Effort  should  continue  to 
understand  gage/propellant  interaction  for  stress  gages  as  well  as  strain 
gages.  Emphasis  on  the  development  of  more  compliant  strain  gages  for 
propellant  is  strongly  recommended.  Innovative  approaches  (such  as  the  use 
of  conductive  polymers  or  gage  structures  fabricated  from  silicon  with  in- 
tegral strain  gages)  should  be  considered  along  with  more  conventional 
approaches  such  as  those  discussed  in  3.4.2. 
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It  is  also  recommended,  because  of  the  excellent  agreement  with  experimental 
results,  that  consideration  be  given  to  the  use  of  the  analysis  techniques 
discussed  in  this  report  for  correcting  strain  gage  data  for  gage/propellant 
interaction  effects.  This  would  provide  a short  term  substitute  for  the 
more  desirable,  long-term  solution  of  developing  "non interacting"  gages. 
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